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1.0  UmODUCTlON 


The  Air  Force  Rocket  Propulsion  Laboratory  (AFRPL)  regularly  conducts 
subscale  thrust  stand  tests  using  the  Ballistic  Test  and  Evaluation  System 
(BATES).  The  tests  serve  to  characterize  propellant  performance  in  standard 
test  motors,  and  the  performance  parameters  usually  measured  are  chamber 
pressure  and  thrust.  Although  these  ballistic  measurements  are  essential 
inputs  in  assessing  propellant  formulations,  they  yield  little  insight  to 
combustion  or  plume  radiation  phenomena. 

Increased  efforts  in  the  development  of  missile  early  warning, 
identification,  and  damage  assessment  systems  require  more  accurately  defined 
plume  radiation  and  particle  property  knowledge.  Plume  signature  measurements 
can  serve  to  help  verify  Joint  Army,  Navy,  NASA  and  Air  Force  (JANNAF)  plume 
signature  prediction  code  output  since  particle  properties  are  sensitive  and 
not  well-characterized  input  parameters  to  these  codes.  Better  quantification 
will  reduce  the  uncertainty  of  code  predictions.  Spectral  signatures  and 
visibility  figures  of  merit  can  also  be  used  to  assess  propellant  formulation 
observability  trade-offs.  The  need  for  plume  measurements  is  especially  clear 
in  i.ight  of  targets  of  opportunity  offered  by  the  ongoing  AFRPL  BATES 
experiments. 

BATES  solid  rocket  motor  exhaust  plume  measurements  were  sought  after 
in  an  in-house  project  entitled  Solid  Rocket  Motor  Test  Analysis.  This  work 
was  directed  _v>ward  three  objectives:  (1)  to  measure  plume  infrared  (IR) 
radiative  emission  and  absorption  characteristics;  (2)  to  measure  plume 
ultraviolet  (UV)  radiation  signatures;  and  (3)  to  measure  plume  particle 
sizes.  An  ancillary  goal  was  to  investigate  use  of  a  measurement  in  serving 
as  an  observability  figure  of  merit,  in  particular,  laser  transmissivity.  The 
original  intent  of  the  measurements  was  to  reduce  particle  sizing  laser 
transmission  and  scattering  measurements  to  particle  size  distribution 
functions  (PSDF’s)  and  to  reduce  IR  emission/absorption  (IR-E/A)  measurements 
to  species  temperature  and  concentration  profiles.  Supplemental  in-vitro 
particle  size  measurements  were  also  sought  after  via  particulate  sampling 
probes.  All  these  measurements  were  to  be  compiled  in  a  data  file  available 
for  predictive  code  and  propellant  formulation  analyses  as  needed.  Actual 


results  fall  somewhat  short  of  these  objectives  and  are  detailed  within  the 
text  of  this  report. 

The  rocket  plume  measurement  instrumentation  was  chosen  to 
specifically  address  specific  measurement  goals.  The  instrumentation  includes 
a  multi-wavelength,  single  line-of-sight  IR-E/A  system,  a  UV  emission 
spectrometer,  an  exhaust  flowfield  particle  collection  probe  and  a  laser 
transmissometer/scattering  system  capable  of  recording  beam  attenuation  and 
Mie  scattered  light  in  both  polarized  states  and  in  its  unpolarized  form.  An 
additional  instrument — a  narrow  field  of  view  radiometer — was  applied  in  this 
project  but  only  in  a  few  firings  late  in  the  experimental  sequence.  This 
application  is  discussed  in  Appendix  D. 

The  specific  data  reduction  and  retrieval  methods  associated  with 
deconvolving  plume  species  temperature  and  concentrations  from  E/A  data  and 
plume  particle  sizes  from  laser  transmission  and  scattering  data  are  described 
in  Sections  3.1  and  3.3,  respectively,  and  in  Appendices  A  and  B.  The  primary 
focus  of  this  report  is  the  AFRPL  rocket  exhaust  measurement  methodology  and 
ensuing  results. 

2.0  APPARATUS 

Descriptions  of  the  facilities  and  instruments  used,  calibration  steps 
applied,  and  the  procedures  for  signal  conditioning  data  output  follows.  The 
reasons  for  applying  this  instrumentation  and  the  associated  signal 
conditioning  and  calibration  techniques  are  discussed  in  this  section. 

2.1  BATES  Test  Stand  and  Motor 

The  plume  measurements  were  conducted  at  an  ambient  test  facility,  Pad 
5a  of  Experimental  Area  1-32.  The  Allegany  Instrument  Company  thrust  stand  at 
Pad  5a  is  designed  to  measure  from  0  to  500  and  from  0  to  10,000  lbs  with 
+0.125%  accuracy.  Data  was  input  digitally  using  the  Nova  DATUM  Data 
Acquisition  and  Control  System  (DACS)  at  a  1000  samples/second  rate  for 
ballistic  motor  parameters  and  at  a  250  samples/second  rate  for  the  optical 
plume  measurements. 


The  15-lb  BATES  motor  used  with  the  plume  measurements  is  particularly 
well-suited  to  tactical  solid  rocket  motors.  Although  the  15-lb  BATES  motor 
uses  different  insulation  and  ignitor,  its  operating  conditions  and  7-in. 
grain  diameter  are  similar  to  those  found  in  a  field  application.  Figure  1  is 
a  schematic  of  a  15-lb  BATES  motor  whose  motor  characteristics  are  summarized 
in  Table  1. 


2.2  Propellants 

The  propellants  used  were  reduced  smoke,  low  burn  rate  formulations, 
and  most  were  targets  of  opportunity  afforded  by  ballistic  evaluations  under¬ 
taken  in  air-launched  propellant  studies  (Refs.  1  and  2).  Aluminum  was  added 
to  the  propellant  more  to  stabilize  combustion  than  as  a  fuel.  Table  2  lists 
ingredient  compositions  of  representative  reduced  smoke  propellants  used  in 
the  firings  reported  here. 

Representative  exhaust  gas  mole  fractions  have  been  calculated  for  a 
16Z  Al  fuel  mass  fraction-loaded  propellant  plume  (designated  Smoky)  (Ref.  4) 
and  for  three  general  reduced  smoke  plumes  classes,  High  Temperature  Nethal 
Fuel  Solid  Rocket  Plumes  (HTMFP:  CO,  HC1,  AI2O3),  Low  Temperature  Nethal  Fuel 
Solid  Rocket  Plumes  (LTMFP:  H2O,  CO,  HC1,  AI2O3),  and  Reduced  Smoke  Low 
Visibility  Solid  Rocket  Plume  (RSP:  H20,  C02,  CO,  HC1,  A1203)  (Ref.  3). 
These  are  listed  in  Table  3.  Note  that  most  plumes  measured  in  this  project, 
including  those  propellants  listed  in  Table  2,  fall  in  the  RSP  class. 


Table  1.  15-Pound  11118  Motor  Configuration 


Nominal  Propellant  Height,  kg 

7 

Maximum  Pressure,  MPa 

14 

Pressure  Range,  MPa 

1  to  14 

Motor  Weight,  kg 

60 

Nozzle  Divergence  Angle,  deg 

15 

Nozzle  Convergence  Angle,  deg 

Grain  Configuration,  cm 

45 

Length 

30 

Outer  Diameter 

17.15 

Inner  Diameter 

11.68 

Heb ,  cm 

2.54 

Web  Action  Time,  s 

1  to  5 

2 

Burning  Surface,  cm 

1320 

Mass  Flow  Rate,  kgs 

1  to  7 

Burning  Rate,  mms  ’l 

3  to  25 

Pressure  Neutrality 

1.08 

Port-to-Throat  Ratio 

3  to  32 

Nozzle  Throat  Diameters,  cm 

Ignitor  BKNO3  pellets 

1  to  6 

Table  2.  Height  Percentage  Composition  of  Representative 
AFRPL  Seduced  Saoke  Propellants 


Ingredient 

ESS 

JAC 

n» 

CHS 

PUT 

Ammonium  Perchlorate  (AP) 

86.0 

86.0 

86.5 

76.0 

75.0 

Hydrocarbon  Binder 

8.64 

10.8 

9.7 

9.7 

9.7 

Dioctylazelate  (DOA) 

2.15 

2.19 

2.0 

2.0 

2.0 

AI2O3 

0.50 

0.50 

— 

— 

— 

Rb  Suppressant 

— 

— 

1.5 

12.0 

5.0 

Dimosyle  Disscynate 

2.19 

— 

— 

— 

— 

Other  Additives 

0.53 

0.53 

0.3 

0.3 

8.3 

The  RSP  formulation  with  its  lov  AI2O3  loading  is  astenable  to  a  laser 
transmissometer/scattering  application  because  multiple  scattering,  focal 


Table  3.  Plume  Species  C<npoaitiou 
of  HTMFP,  LTMFP  and  BSP  Propellants  at  Sea  Level 


Temperature  (K) 

Smoky 

1860 

HTMFP 

1800 

LTMFP 

1000 

ISP 

1000 

Pressure  (atm) 

1.22 

1.0 

1.0 

1.0 

Plume  Model  Fraction 

Gases 

h2o 

0.020 

0.02 

0.30 

0.40 

C02 

0.0037 

— 

— 

0.12 

CO 

0.30 

0.30 

0.02 

0.12 

HC1 

0.13 

0.13 

0.13 

0.18 

h2 

0.41 

0.40 

0.40 

0.40 

n2 

0.068 

0.10 

0.05 

0.03 

Particles 

A1?0-* 

0.07 

0.50 

0.50 

0.05 

(Note  that  H2O,  CO2,  CO  and  HC1 

are  the 

major  UV/IR 

radiating 

RSP  j 

H2  and  N2  are  UV/IR  nonradiating). 


volume  coincidence  and  particle  cross-sensitivity  effects  are  mitigated  by  the 
relatively  low  optical  thickness  of  these  propellant  plumes.  In  contrast, 
highly  aluminized  plumes  from  high  performance  propellants  (HTMFP  and  LTMFP) 
attenuate  the  lasex  signal  to  the  point  where  background  radiation,  multiple 
scattering,  coincidence  and  cross-sensitivity  error  sources  dominate.  These 
issues  are  addressed  in  Section  2.3  and  in  Appendix  A. 


2.3  Particle  Sising  Instrwentation 


Both  nonintrusive  and  intrusive  particle  sizing  instruments  were 
used  in  this  series  of  BATES  motor  plume  measurements.  The  nonintrusive 
instruments  were  optical  and  consisted  of  a  laser  transmission/scattering 
experiment.  Intrusive  techniques  included  probe  sampling  of  the  plume  exhaust 
flow. 


•  V  V  . 
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2.3.1  Laser  Transmission/Scattering  Instruments 


Description  -  The  AFRPL  laser  transmission/scattering  configuration 
was  modeled  after  an  Arnold  Engineering  Development  Center  (AEDC)  particle 
size  retrieval  experiment  with  a  short-duration  (130  ms)  arc-heated  flowfield 
(Ref.  5).  The  intent  was  to  translate  the  AEDC  particle  size  measurement 
method  to  AFRPL  measurements.  This  technique  had  never  before  been  applied  to 
rocket  exhaust  flowfields. 

The  laser  transmission/scattering  measurement  system  orientation  is 
depicted  in  Figures  2  and  3.  The  system  was  designed  to  accomplish  plume 
measurements  in  the  relatively  short-lived  and  hostile  BATES  motor  test 
environment.  The  anticipated  severity  of  vibration  and  temperature  dictated 
installing  a  massive  250-kg  support  structure  to  enclose  the  measurement 
instruments  and  to  maintain  alignment  between  the  laser  and  detectors.  The 
structure  was  fabricated  from  30-cm  channel  iron  and  was  mounted  in  a  concrete 
foundation.  A  laser  and  detectors  were  attached  to  this  rigid  structure  (Fig. 
3).  A  3.5-cm  thick  rigid  aluminum-backed  foam  insulated  the  entire  structure 
and  provided  the  enclosed  instrumentation  with  some  measure  of  protection  from 
both  the  hostile  exhaust  plume  environment  and  shrapnel  from  occasional  nozzle 
failures.  The  structure  was  set  2  meters  axially  from  the  exit  plane  of  a 
mounted  15-lb  BATES  motor  casing.  The  entire  system  was  kept  indoors  to 
maintain  a  controlled  ambient  temperature. 

Figure  2  is  a  schematic  of  the  instrumentation  system  which  starts 
from  a  4  Watt  Argon-ion  laser  located  across  the  plume  from  a  laser  power 
meter;  this  is  the  laser  transmissometer  portion  of  the  experiment.  There  are 
six  detectors  arranged  in  a  plane  perpendicular  to  the  plume  axis  and  a 
seventh  broadband  detector  for  scattering  measurements.  Electro-optical 
parameters  associated  with  this  system  (Table  4)  are  detailed  in  the  ensuing 
discussion. 

The  Control  Laser  Corporation  Model  552A  4W  Argon-ion  laser  operates 
multi-line,  but  delivered  1.5  W  average  power  at  514.5  nm  for  these 
measurements.  The  beam  was  interrupted  for  approximately  0.2  s  midway  through 


each  motor  burn  with  a  solenoid-operated  shutter  to  enable  subsequent 
correction  of  power  meter  and  scattering  detector  output.  The  output  at  the 
laser  source  was  independently  monitored  to  spot  and  correct  laser  power 
fluctuations  and  long  term  drift. 

During  operation  the  laser  light  first  traversed  a  45°  polarizer,  thus 
simulating  unpolarized  light  (comprised  of  equal  parallel-  and  perpendicular- 
polarized  components).  The  United  Detector  Technology  (UDT)  Model  21A  optical 
power  meter  used  to  detect  the  transmitted  laser  signal  had  a  broadband  (380- 


Figure  3.  AFRPL  Laser  Transmission/Scattering 
Measurement  System  During  Motor  Firing 

1100  nm)  spectral  range  silicon  PIN  photovoltaic  photodiode  detector.  A  10-nm 
full  width  half  maximum  (FWHM)  bandwidth  filter  centered  at  the  514.5-nm  laser 
wavelength  was  placed  between  the  power  meter  detector  and  a  UDT  Model  2500 
integrating  sphere  exit  port  and  enabled  full  use  of  the  power  meter's  0.1  pW 
-  10W  sensitivity  range. 

The  six  scattering  detector  angles  are  indicated  on  Figure  2  as  15, 
30,  50,  90,  120,  and  150  degrees  from  the  forward  scattering  direction.  The 
scattering  detection  assemblies  at  each  of  these  angles  were  identical,  each 
with  10-cm  diameter  f/2.0  lenses.  Figure  4  shows  the  configuration  of  the 
laser  scattering  detector  and  Figure  5  is  a  photograph  of  the  laser  scattering 
detector  cross  section. 

These  scattering  detectors  were  designed  to  operate  in  the  same  plane 
orthogonal  to  the  plume  axis.  All  six  were  to  be  equidistant  to  the 
centerline,  defining,  along  with  the  laser  beam  shape  profile,  a  scattering 
volume  on  which  ensemble  Mie  particle  deconvolution  techniques  could  be 
applied.  Application  of  the  Mie  theory  and  the  assumptions  used  are  described 
in  Appendix  A.  AFRPL  data  application  retrieval  methods  are  discussed  in 
Section  3.1. 
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Table  4 


.  Soae  Laser  Tranaaias ion/Scattering 
Electro-Optical  Paraaeters 

i-Asm 

Beam  Power  (TEMqq)  ■  4W  nominal,  1.5  W  @  514.5  nm 
Beam  Width  ■  2.0  mm 
Beam  Shape  *  Gaussian 
Wavelength  ”  514.5  nm 

POWER  METER  DETECTOR  (Si  Photodiode) 

Optical  Bandwidth  *  10  nm  centered  @  514.5  nm 

Electrical  Bandwidth  *  DC  to  500  kHz  (but  with  total  power  meter 
system  at  DC  to  200  Hz) 

2 

Detector  Active  Area  **  1.0  cm 

Noise  Equivalent  Power  (NEP)  ■  5.9  *  10  ^  W  (1000  Hz,  1  Hz,  880  nm) 

SCATTERING  DETECTORS  (Si  Photodiode) 

Optical  Bandwidth  ■  1  nm  centered  @  514.5  nm 

Electrical  Bandwidth  *  DC  to  0.1  MHz 

2 

Detector  Active  Area  ■  1.0  cm 

NEP  -  1  x  10-12  W  (1000  Hz,  1  Hz,  850  nm) 


Each  of  the  six  scattering  detectors  had  four  detection  modes: 
vertically  polarized,  horizontally  polarized,  nonpolarized  and  dark.  For  each 
detector,  rapid  sequencing  of  18  polarizing  filters  was  initially  used  to  give 
100-ms  time  resolution  for  each  of  the  detection  modes,  but  because  of 
stresses  induced  by  this  dynamic  mechanical  environment,  the  improvements 
described  by  Reference  6  were  implemented.  With  these  improvements  time 
resolution  still  remained  at  100  ms,  but  12  aluminum-reinforced  DC  solenoid- 
actuated  detectors  sequencing  every  400  ms  now  comprised  the  polarization 
system.  Figure  6  illustrates  how  a  scattering  detector  channel  signal  with 
such  a  sequence  was  designed  to  appear.  There  was  an  approximately  10-ms 
settling  time  in  the  filter  activation  sequence.  This  sequence  was  similar  to 
the  pre improvement  filter  sequence,  with  differences  accounted  for  by  data 
preprocessing. 

The  six  scattering  detectors  were  UDT  Model  500D  silicon  photovoltaic 

-3  -12 

detectors  responding  from  300  to  1100  nm  with  10  to  10  W  selectable 

nominal  sensitivity  at  850  nm.  A  narrow  0.1-nm  FWHM  interference  filter 
separated  the  scattered  laser  light  from  other  plume  emission  wavelengths. 
The  seventh  detector  was  also  a  UDT  Model  500  photodiode,  but  it  had  a 
relatively  broadband  filter  (10-nm  FWHM  centered  at  514.5  nm).  It  was  mounted 
alongside  the  90°  scattering  angle  detector  and  placed  in  an  unfocused 
detector  assembly  directed  upstream  from  the  laser  beam  but  still  at  the  plume 
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centerline.  It  wee  to  serve  as  a  plume  radiation  correction  signal  for  the 
six  narrowband  scattering  detector  outputs. 


DCTECTOH 
SIGNAL 
(m  volts) 


TIME  |m  sec) 

Figure  6.  Idealized  Scattered  Light  Data  Polarizer 
Activation  Sequence 


Signal  Conditioning  -  Signals  from  the  laser  power  meter  and  seven 
scattering  detectors  were  amplified  at  each  instrument  with  a  preamplifier 
then  carried  on  300  ft  of  coaxial  shielded  cable  to  the  data  acquisition 
system.  Signals  of  0-10  mv  amplitude  were  digitized,  then  recorded  at  250 
samples/s  on  a  100  Hz-filtered  NOVA  DATUM  digital  data  acquisition  system 
magnetic  storage  tape. 

Calibration 


A.  Laser  and  Laser  Power  Meter  -  The  laser  power  meter  operated  in  an 
uncalibrated  absolute  responsivity  mode  since  transmissometer  data  reduction 
was  accomplished  mainly  by  measuring  laser  transmittance  values  during  rocket 
motor  firings.  Since  the  zero  transmittance  level  was  accomplished  at 
midfiring  by  a  laser  beam  breaker,  background  levels  were  accounted  for  and 
did  not  need  to  be  calibrated  out.  However,  two  important  calibration 
considerations  still  remained:  verification  of  laser  power  stability  and 
background  and  random  noise  determination. 


A  data  channel  from  the  laser  source  was  directly  monitored  to  give 
continuous  laser  output  voltage  vs  time.  This  data  was  divided  into  two  equal 
time  intervals  and  histograms  depicting  the  laser  output  voltages  for  each  of 
these  intervals  were  calculated,  then  overlaid  with  theoretical  normal 
distribution  curves  determined  from  data  means  and  standard  distributions.  If 
the  laser  noise  was  indeed  Gaussian  random  in  form  (and,  therefore, 
subtractable  as  white  noise),  the  histogram  bins  would  follow  the  normal 
distribution  model;  moreover,  there  would  be  little  variation  between  the  two 
independent  time  intervals.  An  application  of  this  analysis  to  a  BATES  firing 
is  depicted  by  Figures  7-9. 


I  CU  DARK  CURRENT  TEST  23  MAR  84 
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Figure  7.  Laser  Output  Voltage  vs  Time 


Note  that  the  quantization  observed  in  Figure  7  is  a  consequence  of 
the  low  signal  rate  from  the  laser  output  voltage  channel.  The  data  steps 
recorded  in  the  figure  are  highly  resolved  voltage  counts  from  the  data 


PROBABILITY  DENSITY 


acquisition  system.  This  quantization,  when  coupled  with  the  small  standard 
bin  widths  used,  explains  the  gaps  in  the  histograms  shown  in  Figures  8  and  9. 


MEAN  SQUARE 
460.633 


RMS  VARIANCE  CHI  SO 
21.923  5*996  843*92 


HERN  VALUE 
21*766 


STO  DEV 
2*449 


hlNIlUlh  H  AX  I  HUH  VALID  OBS 

12*733  27.823  1358  /V 


UNFILTERED  ORTA 


Figure  8.  Probability  Density  vs  Laser  Output  Voltage,  Interval  1 


To  determine  the  noise  contribution  from  any  coherent  components,  a 
power  spectrum  analysis  was  also  undertaken.  Figures  10  and  11  show  smoothed 
Hamming  window  power  spectra  of  the  laser  voltage  output  data  given  in  Figure 
7.  These  were  taken  from  the  same  two  time  intervals  the  histograms  of 
Figures  8  and  9  were  drawn  from. 


For  the  firings  reported  herein,  the  laser  output  voltage 


noise  was 


PROBABILITY  DENSITY 
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MEAN  VALUE 
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STO  DEV 
1 .797 
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Figure  9.  Probability  Density  vs  Laser  Output  Voltage,  Interval  2 

found  to  be  incoherent,  making  the  laser  transmissometer  data  suitable  for  the 
time-averaged  data  analysis  procedure  explained  in  Section  2.4. 

The  analysis  of  background  and  random  noise  levels  of  the  laser  power 
meter  was  accomplished  in  a  similar  manner.  The  laser  power  meter  output 
voltage  channel  was  periodically  monitored  for  random  noise,  coherent  noise 
forms  and  responses  to  leakages  or  internal  radiation.  This  monitoring  was 
conducted  via  a  dark  current  measurement  calibration  step  where  the  power 
meter's  integrating  sphere  input  aperture  was  blocked  with  a  light-tight 


BATES  2628-005  28/APR/82  TIME  INTERVAL  -.4  TO  2.8  SECONDS  FROM  IGNITION 
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Figure  10.  Laser  Output  Voltage  Power  Spectrum,  Interval  1 

cover.  As  with  the  previous  analysis,  the  next  goal  was  to  identify  noise 
forms  present  in  the  calibration  data  and  to  eliminate  colored  noise 
constituents.  Results  of  a  dark  current  calibration  and  histogram  analysis 
are  shown  in  Figures  9  and  10.  Power  spectra  are  given  in  Figure  11. 


Note  the  quantization  effect  evidenced  in  Figure  12.  As  with  Figure  7 
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Figure  11.  Laser  Output  Voltage  Power  Spectrum»  Interval  2 

this  is  a  result  of  the  disparity  between  the  low  signal  rate  and  high  data 
acquisition  system  resolution.  Empty  histogram  bins  of  Figures  13  and  14 
again  result  from  this  quantization. 

Note  that  60  and  68  Hz  noise  spikes  stand  out  in  the  otherwise 
relatively  incoherent  noise  in  the  laser  power  meter  dark  current  calibration 


Figure  12.  Laser  Power  Meter  Output  Voltage  vs  Time 


cut  shown  in  Figure  15.  The  60  Hz  noise  has  some  background  AC  current 
crosstalk  and  its  elimination  requires  searching  characteristic  AC  sources. 
The  68  Hz  spike,  however,  is  more  difficult  to  control  and  a  narrow  band  notch 
filter  may  be  the  best  fix. 


®*  Scattering  Detectors  -  Calibrating  the  absolute  responsivity  of 
the  scattering  detectors  presented  a  difficult  problem.  The  detectors  were 
essentially  locked  into  the  scattering  geometry  depicted  in  Figure  2,  and  it 
was  thought  that  removing  them  from  the  laser  scattering  system  for  absolute 
calibration  with  a  blackbody  reference  source  would  prohibitively  alter  the 
firing-to-firing  alignment.  However,  minor  adjustments,  especially  along  the 
plume  radius,  were  relatively  easy  to  do.  To  maintain  the  alignment 
integrity,  a  quartz-halogen  lamp  was  placed  at  the  focal  volume  defined  by  the 
scattering  detectors.  The  lamp  filled  the  detectors'  FOVs.  This  absolute 
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Figure  13.  Probability  Density  vs  Laser  Power, 
Meter  Output  Voltage 


calibration  technique  assumes  spatiel  radiance  uniformity  across  the  lamp, 
equal  scattering  detector  instantaneous  fields  of  view  (IFOVs)  and  tbe  same 
relative  angle  between  the  lamp  and  each  detector. 

To  ensure  the  integrity  of  the  responsivity  calibration  several 
supporting  calibration  parameters  were  necessarily  determined  to  fully 
characterize  influencing  error  source!*  (Ref.  7).  These  are: 


(1)  absolute  spectral  positioning; 

(2)  spectral  resolution; 

(3)  out-of-band  rejection; 

(4)  frequency  response  determination; 


SCALAR  VALUE 

Figure  14.  Probability  Density  vs  Laser  Power  Meter 
Output  Voltage 


(5)  FOV  mapping;  and 

(6)  background  and  noise  level  determination. 

These  factors,  and  especially  Factors  5  and  6,  are  discussed  in  further  detail 
in  this  section. 


Out-of-band  rejection  was  also  easily  accounted  for  since  the  Si 
photodiode  response  cut  off  above  1.06  pm,  and  since  solid  rocket  motor  plume 
gaseous  species  radiate  virtually  negligible  amounts  within  the  Si  spectral 
response  (Ref.  8).  Most  radiation  at  these  optical  frequencies  occurs  with 
hot  plume  particle  continuum  radiation  in  addition  to  particle-scattered  and 
particle-transmitted  light  from  background  sources.  Figure  17  shows  that  even 
with  2000K  particles  emitting  as  a  blackbody  continuum,  there  is  not  much 
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Figure  15.  Laser  Power  Meter  Output  Voltage 
Power  Spectrum,  Interval  1 


radiative  flux  in  the  vicinity  of  0.5  Mm.  This  is  further  substantiated  by 
the  Naval  Postgraduate  School's  experience  with  dif fractively  scattered  light 
measurements  (Ref.  9).  Hence,  out-of-band  leakage  is  not  a  problem  for  the 
scattering  detector — or  for  that  matter,  laser  power  meter — application. 
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Figure  16.  Laser  Power  Meter  Output  Voltage 
Power  Spectrum,  Interval  2 

The  Si  photodiode  scattering  detector  10-90%  rise  time  was  2  x  10~7s. 
Further  scattering  detector  electro-optical  parameters  were  given  earlier  in 
Table  4. 


This  brings  us  to  the  FOV  mapping  step  in  the  scattering  detector 


Wavelength  (|im) 


Figure  17.  Blackbody  Spectral  Radiant  Exitance  Distribution 
Curves  for  Temperatures  Between  1000  and  2000K 

calibration  sequence.  Ideally,  a  polar  plot  of  the  two-dimensional  scattering 
detector  radiometer  FOV  would  be  performed  to  quantify  the  instrument  response 
as  a  function  of  incidence  angle.  This  is  essential  to  ascertain  that  no 
detector  spatial  detectivity  anomalies  exist  and  to  define  the  normalization 
curves  which  characterize  the  spatial  radiation  distribution  within  each 
scattering  detector's  focal  plane.  This  polar  plot  should  be  conducted  with  a 
collimated  beam  of  sufficient  intensity  and  small  enough  area  to  elicit  a 
point-by-point  detector  system  response. 

The  facilities  to  accomplish  this  (e.  g.,  a  long  focal  length 

» 

collimator  system)  were  unavailable  to  AFRPL  at  the  time  the  BATES  motor 
measurements  described  in  this  report  were  taken.  In  lieu  of  the  polar  plot, 
a  one-dimensional  FOV  mapping  calibration  along  the  plume  radial  direction  was 
performed  for  each  firing.  Because  the  detectors  were  intended  to  measure 
angular  scattering  intensities  in  the  radial  direction,  the  second  dimension, 
plume  axial  direction,  was  determined  to  be  of  lesser  importance.  Detector 
response  was  measured  as  a  function  of  distance  along  the  plume  radius.  The 

mapping  served  as  an  alignment  and  focusing  aid,  but  the  scattering  detector 

assembly  afforded  limited  adjustment  in  the  plume's  radial  direction. 

The  background  and  noise  determination  step  in  the  calibration 
sequence  was  identical  to  that  just  reported  in  the  laser  power  meter 


calibration.  Results  o£  one  such  scattering  detector  dark  current  calibration 
are  given  in  Figures  18-22.  As  with  Figures  7  and  12,  quantization  effects 
are  shown  to  be  present  due  to  a  difference  between  the  signal  rate  and  data 
acquisition  system  time  resolution. 


CU  DARK  CURRENT  TEST  23  MflR  84 
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Figure  19.  Probability  Density  vs  Scattering  Detector  5 
Output  Voltage  vs  Tine,  Interval  1  V 

% 

and  is,  to  collect  unbiased  solid  rocket  motor  exhaust  particles  to  verify 
motor  performance  and  plume  signature  prediction  codes.  The  same  aim  was  an 
ancillary  objective  of  the  present  AFRPL  effort,  but  a  further  goal — to 
correlate  nonintrusive  laser  particle  size  measurements  with  probe-captured 
particle  data — was  sought.  The  now  more  pressing  need  to  collect  unbiased  and 
undisturbed  particles  both  in  terms  of  size  and  surface  morphology  to  verify 
the  code  gave  emphasis  to  designing  a  supersonic  shock-capturing  probe.  It  is 
hoped  that  the  new  design  will  capture  unbiased  rocket  exhaust  particles  in 
the  field  which  has  never  been  accomplished  to  any  significant  degree. 
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Figure  20.  Probability  Density  vs  Scattering  Detector  5 
Output  Voltage  vs  Time,  Interval  2 


At  least  four  basic  problems  oust  be  solved  to  successfully  extract  an 
unbiased  sample  of  particulates  from  a  rocket  exhaust.  They  are: 


(1)  the  collection  probe  may  bias  the  particle  diameter  measurement 
via  bow  shock  and  recirculation  effects,  changing  the  collected  particle  size 
distribution; 


(2)  particulate  entrainment  problems  may  introduce  foreign  particles 
into  the  flow; 


(3)  agglomerates  of  smaller  particles  may  be  broken  up  during 
collection  and  handling;  and 
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Figure  22.  Scattering  Detector  5  Output  Voltage 
Power  Spectrum.  Interval  2 


A  major  advantage  of  intrusive  collection  techniques  over  the  nonintrusive 
optical  particle  sizing  means  is  that  heavily-laden  particulate  flows  can  be 
investigated;  optical  thickness  is  of  no  concern. 


Two  intrusive  flowfield  particle  collection  probes  are  described  in 
this  section.  The  first,  a  subsonic  particle  collection  probe,  has  been 
applied  to  a  number  of  the  firings  described  in  this  report.  The  latter,  a 


27 


shockless  supersonic  inlet  particle  collection  probe,  has  gone  through  its 
final  design  phase  under  the  AFRPL  particle  sizing  effort  reported  here. 

Subsonic  Particle  Collection  Probe  Description  -  The  AFRPL  subsonic 

2 

particle  collection  probe  is  the  36.5  x  5  cm  torpedo-shaped  aspirated  tube 
shown  in  Figure  23.  In  its  current  application  the  probe  is  placed  on  the 
plume  centerline  and  approximately  seven  meters  downstream  from  the  BATES 
nozzle  exit  plane.  A  portion  of  the  subsonic  plume  flow  is  captured  within 
the  0.5-cm  entrance  orifice  and  particles  within  this  flow  adhere  to  a  number 
of  25-micron  platinum  wires  strung  in  an  assembly  within  the  probe.  Collected 
particles  were  then  submitted  for  further  analysis  to  a  Nanometrics  Model  HPS 
70  Field  Emission  Scanning  Electron  Microscope  capable  of  70&  resolution 
photomicrographs.  A  Bausch  and  Lomb  QMS  image  analyzer  was  also  applied  but 
with  limited  success,  as  described  later  in  Section  3.1. 


Figure  23.  Subsonic  Particle  Collection  Probe 

The  platinum  wire  collection  assembly  shown  in  Figure  24  represents  an 
improvement  implemented  at  AFRPL  in  July  1982.  Prior  to  this,  particles 
■within  the  probe  impacted  on  a  blunt  2-cm  diameter  aluminum  plug.  When  placed 
in  the  plume  flow,  a  stagnation  point  was  created  at  the  center  of  the  plug, 

* 

.with  most  of  the  particle-laden  flow  moving  axisymraetrically  around  the 
•  collection  surface,  biasing  the  collection  away  from  the  smaller  particles 
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Figure  24.  Platinum  Wire  Collection  Assembly 


which  tend  to  follow  the  flow  (Ref.  12).  The  small  cross  sectional  area  of 
the  platinum  wires  alleviates  much  of  this  bias. 

Despite  the  improvement,  collected  particles  still  tend  to 
inaccurately  represent  plume  flowfield  particles  because  of  collection 
methodology  error  sources  such  as  particle  agglomeration,  bow  shock,  platinum 
wire  impact,  and  particle  handling  effects.  This  points  to  the  difficulty  in 
accurately  applying  a  sampling  technique  to  rocket  exhaust  field  measurements. 


Supersonic  Shock-Swa  I  1  >iwi  ng  (Jucuth  Probe  Descri.pl  ion  -  Although  it  is 
virtually  impossible  to  eliminate  all  ot  tin  effeits,  we  ran  minimize  them  and 
construct  a  probe  so  that  the  extent  ut  bias  can  be  estimated. 


The  supersonic  particle  collection  device  has  three  objectives: 

(1)  to  collect  contaminated  particles  for  chemical  and  surface 
property  determination; 

(2)  to  collect  a  particle  sample  whose  total  mass  is  not  biased  by 
the  collection  scheme;  and 

(3)  to  collect  a  particulate  sample  where  particles  do  not  pass 
through  a  normal  shock  which  will  help  determine  if  particle  agglomeration 
occurs  to  a  statistically  significant  extent. 

The  particle  capture  probe  shown  schematically  in  Figure  25  attempts 
to  address  all  these  goals.  The  probe  consists  of  a  sharp  tungsten  tip 
designed  to  efficiently  capture  or  swallow  a  stream  tube  of  one  gram  per 
second  of  plume  in  a  supersonic  flow  if  the  internal  pressure  is  kept  low. 
Similar  tips  designed  by  one  of  the  authors  at  AEDC  have  been  exposed  to  the 
supersonic  flow  of  arc  heated  air  at  envelopes  as  high  as  2000  B/lb  and  pitot 
pressures  as  high  as  80  atm  for  periods  of  1.0  second  with  no  apparent  ill 
effect.  The  performance  of  the  tungsten  tip  in  the  highly  erosive  rocket 
engine  particulate  flow  is  not  expected  to  be  as  good  as  has  been  obtained  in 
"clear"  air,  so  the  tip  has  been  made  removable  to  facilitate  its  replacement 
and  refurbishment. 

To  ensure  that  the  probe  does  not  stagnate  the  captured  stream  tube,  a 
stream  of  cold  nitrogen  has  been  introduced  to  act  as  an  ejector  flow,  thus 
reducing  the  tip  exit  pressure  to  7.0  psia.  The  probe  tip  will  be  pumped 
through  when  the  nitrogen  is  turned  on.  The  plume  flow  ratio  of  14:1  was 
selected  to  decelerate  and  quench  the  hot  rocket  motor  plume  particles.  The 
specimen  temperature  is  reduced  to  about  270  R  to  preclude  burning  the  filter. 
The  mixing  passage  flow  rates  and  areas  were  sized  to  prevent  formation  of  a 
normal  shock  at  the  design  conditions.  A  series  of  oblique  shocks  will 
decelerate  this  combined  flow  to  a  slightly  lower  Mach  number.  A  normal  shock 
is  expected  to  raise  the  static  pressure  to  about  the  30  psia  required  to 
force  the  captured  rocket  plume  particles  through  the  0.025  ym  removable 


0-025  MICRON  FILTER 


millipore  filter.  A  pressure  passage  has  been  provided  to  measure  the  gas 
static  pressure  at  the  point  that  the  core  flow  discharges  into  the  cold 
quench  stream. 


The  removable  raicropore  filter  has  been  selected  to  collect  up  to  a 
gram  of  particulate  without  creating  undue  filter  pressure  loss.  The  filter 
will  be  analyzed  to  determine  the  capture  particulate  diameter  distribution. 
Since  the  expansion  field  will  tend  to  separate  the  particulate  as  a  function 
of  diameter,  a  careful  statistical  analysis  of  the  filter's  PSDF  will  be  made 
to  determine  the  effect  of  the  expansion  process  on  the  particulate  input 
location. 

Since  the  probe  is  designed  to  stay  in  the  rocket  plume  flow  for  only 
1.0  second,  we  have  also  designed  a  probe  injection  mechanism  and  shield.  In 
practice  the  injector  will  rotate  the  probe  into  the  plume  flow  during  the 
motor  firing  to  avoid  the  debris  and  ignitor  smoke  from  the  ignition 
transient.  A  shield  will  be  rotated  into  the  flow  along  with  the  probe  to 
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ensure  it  is  only  exposed  to  the  plume  flow  at  a  predetermined  location  in  the 
plume.  A  pair  of  pneumatic  actuators  has  been  configured  to  rotate  the  probe 
and  shield  into  and  out  of  the  flow  and  to  cover  and  uncover  the  particle 
collection  probe  in  the  desired  sequence.  The  probe  stand  is  designed  to  be 
adjustable  to  accommodate  a  variety  of  rocket  motor  test  stands  at  AFRPL. 

The  particle  collection  test  sequence  will  begin  by  first  collecting 
particles  in  the  least  severe  area  of  the  plume;  thus  our  first  objective  of 
collecting  uncontaminated  particles  will  be  achieved  with  little  risk.  On 
subsequent  tests  the  probe  will  be  located  upstream  of  test  nozzle  Mach  disks 
to  obtain  a  sample  with  subsonic  plume  flow.  If  sufficient  progress  has  been 
made,  the  probe  will  be  redesigned  internally  to  produce  a  completely 
isentropic  internal  flowfield  to  decelerate  the  particulates  without  a  normal 
shock.  This  is  so  particulate  agglomeration  can  be  studied. 

In  addition  to  collecting  particles  for  later  analysis,  the  particle 
probe  can  be  used  to  determine  the  particle  mean  flux  and  PSDF  maps  of  the 
entire  plume  which  can  then  be  checked  against  the  flowfield  predictions.  We 
can  also  obtain  data  on  the  distribution  of  impurities  injected  into  the  flow 
from  ablative  nozzle  walls.  The  determination  of  the  impurities  such  as 
carbon  particles  could  be  expected  to  have  a  large  influence  on  plume 
radiation  and  may  explain  what  has  heretofore  been  consigned  as  anomolous 
plume  radiation  data. 

2.4  Plume  Signature  Instrumentation 

Plume  signature  instrumentation  consisted  of  infrared  and  ultraviolet 
scanning  monochromators  used  in  application-specific  configurations.  These 
will  be  described  in  this  section. 

2.4.1  Infrared  Emission/Absorption  (IR-E/A)  Instrumentation  -  tR-E/A 
measurement  techniques  have  been  successfully  applied  to  yield  specie 
temperatures  and  partial  pressure  profiles  with  long  duration  single-phase 
liquid  rocket  engines  (Refs.  13  and  14),  but  the  application  to  short  duration 
two-phase  solid  rocket  motors  is  unique  to  this  AFRPL  experiment.  Because  of 
the  short  (2-5  s)  15-pound  BATES  firing  times,  it  was  found  necessary  to 
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abandon  the  traditional  IR-E/A  multiposition  spatial  scanning  technique 
described  in  Reference  13  and  to  work  from  a  fixed  position  at  the  plume 
centerline  near  the  motor  exit  plane.  To  compensate,  it  was  decided  to  obtain 
multispectral  data  at  the  fixed  IR-E/A  station. 

The  physics  behind  the  IR-E/A  temperature  and  partial  pressure 
inversion  scheme  is  introduced  in  Appendix  B  of  this  report.  The  inversion 
technique  adopted  by  AFRPL  is  described  in  Section  3.3  of  this  report. 

Description  -  Schematics  of  the  system  layout  with  respect  to  the 
BATES  motor  and  plume  are  shown  in  Figures  26  and  27  with  Figure  26  showing 
the  top  view  relative  to  a  BATES  15-pound  motor  and  Figure  27  representing  the 
aft  end  perspective. 


^  FAST  SCAN 
MONOCHROMATOR 


Figure  26.  Top  View  of  IR-E/A  System 


The  radiation  source  system  consited  of  a  1-in  aperture  Barnes  Model 
11-210  blackbody  radiation  source  set  at  1000C.  The  source  radiation  was 
chopped  at  1670  Hz  and  was  focused  to  the  plume  centerline  by  a  10-cm  focal 
length  f/2.5  fused  silica  piano  convex  lens.  The  system  is  shown  in  Figure 
27.  A  housing  which  could  be  flushed  with  gaseous  nitrogen  protected  the 
system. 


The  receiving  system  had  a  Spex  Minimate  Model  1670  grating 
monochromator  with  a  0.22-m  focal  length  and  f/4.0  aperture  ratio,  configured 
in  a  Czerzy-Turner  layout.  The  monochromator  was  equipped  with  a  300 
grooves/mm  diffraction  grating  blazed  to  operate  with  peak  efficiency  at  the  2 
pm  wavelength.  As  was  the  source  system,  a  PAR  Model  125A  optical  chopper  was 
in  the  optical  train  operating  at  670  Hz.  Condensing  optics  consisted  of  the 
identical  fused  silica  lens  found  with  the  source  system.  Again,  a  protective 
housing  with  a  fused  silica  window  surrounded  the  system. 


The  monochromator  was  equipped  with  18  mm  high  and  1.25  mm  wide 

o 

entrance  and  exit  slits  to  provide  50  A  resolution  radiance  to  an  ellipsoidal 
mirror  attachment.  The  mirror  focused  the  exit  slit  image  onto  a  Santa 
Barbara  Research  Center  Model  70750  C126  InSb  detector.  This  detector  later 
replaced  the  original  Spex  Model  1428  thermopile  detector  mainly  because  its 
sub-  ys  rise  time  (vs  the  thermopile's  maximum  4-ms  rise  time)  enabled  it  to 
adequately  resolve  the  signals  received  from  the  fast  grating  scan. 

The  monochromator  was  modified  to  scan  at  9  ym/s  over  its  design  scan 
speed  of  0.013  ym/s.  To  accomplish  this,  a  cam  was  placed  over  the  existing 
lead  screw  and  attached  to  a  sine  bar  which  drove  the  kinematic  grating  mount 


quicker  than  originally  designed.  The  linear  motion  was  preserved  in  the  new 
design.  A  linear  potentiometer  was  fitted  to  the  drive  mechanism  to  produce  a 
voltage  signal  proportional  to  wavelength.  A  24  V  DC  motor  was  also 
attached  to  the  drive  mechanism  to  permit  the  fast  scan. 

o 

The  receiver  system  was  capable  of  a  50  A  resolution  scan  from  1800  to 
o  o 

4900  A  in  350  ms  but  operated  from  2750  to  4250  A. 

Signal  Conditioning  -  The  signal  conditioning  system  is  diagrammed  in 
Figure  28.  The  preamplifier  outputs  an  AC  signal  proportional  to  the 
radiation  impinging  the  InSb  detector.  This  detector  signal  was  then  carried 
through  300  ft  of  coaxial  cable  and,  along  with  the  source  and  receiver 
chopper  synch  channels,  fed  through  two  PAR  Model  128A  lock-in  amplifiers. 
One  lock-in  amplifier  was  synchronized  by  a  1603  Hz  source  chopper  and  the 
other  by  the  630  Hz  receiver  chopper.  In  this  way  the  analog  transmitted 
source  signal  Vf  and  the  radiance  signal  Vg  were  separated  and  then  recorded 
simultaneously.  These  signals,  along  with  the  grating  drives  potentiometer 
signal  ,  were  recorded  on  magnetic  tape  at  250  samples/s  using  a  NOVA  Datum 
digital  data  acquisition  system. 


i  ^  too  ™  ^  DIGIT  It 
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Calibration  -  Absolute  spectral  positioning  was  initially  accomplished 

by  inserting  a  1200  grating/nn  ultraviolet/visible  grating  in  place  of  the  300 

grating/mm  grating  used  in  the  measurement.  A  low  pressure  mercury  lamp 

source  was  then  focused  into  the  monochromator's  entrance  slit  and  the 

motorized  scan  drive  was  actuated.  The  voltage  applied  to  the  drive 

mechanism's  potentiometer  V£n  was  recorded  along  with  the  potentiometer  wiper 

voltage  V.  and  the  detector  response  voltage  Vg.  Characteristic  Hg  line  peaks 
A  o 

at  2536.52,  3131.55,  3663.28  and  4046.56  A  were  then  discriminated  in  a  Vg(V^) 
plot.  Since  the  diffraction  grating  ruling  ratio  of  4  holds  with  a  1200 
gratings/mm  substitution  for  the  300  gratings/mm  calibrating  grating,  it 
follows  that  Vg  peaks  occur  at  V^'s  corresponding  to  1.01,  1.25,  1.47  and  1.62 
pm.  Although  these  wavelengths  are  at  the  low  end  of  both  the  1200 

gratings/mm  grating  efficiency  curves  and  the  InSb  spectral  response,  they 
sufficed  as  initial  curve  fit  points  for  the  function. 

A  “  a  (Vx/VR)  -  b  (1) 

The  wavelength  A  ,  determined  by  the  particular  calibration  a  and  b  curve  fit 
constants  and  output  voltages  and  Vg,  gave  the  absolute  spectral 

positioning  for  the  corresponding  BATES  firing  data.  Implicit  in  this  linear 

form  curve  fit  is  the  assumption  of  a  linear  relationship  between  V  and  Vg 

A 

(i.  e.,  the  potentiometer's  wipe  voltage  linearly  corresponds  to  likewise 
linear  drive  mechanism  and  kinematic  grating  mount  motions). 

To  eliminate  the  effects  of  radiation  from  secondary  orders  of 
diffraction,  three  long  pass  blocking  filters  set  at  1.5,  2  and  3.5  pm  were 
designated  to  be  set  in  the  monochromator's  optical  path  with  cut-on 
wavelengths  depending  on  the  desired  scan  spectral  interval.  Only  radiation 
from  the  first  order  free  spectral  range  was  allowed  through. 

No  F0V  mapping  was  undertaken  with  this  instrument  because  the 
necessary  collimator  was  unavailable  to  AFRPL  at  the  time  these  measurements 
were  taken.  In  addition,  spatial  modulation  of  the  signal  was  assumed  to  be  a 
relatively  insignificant  error  source  because  of  the  characteristic  uniform 
response  over  the  InSb  detector's  small  active  detection  surface  area. 


As  noted  previously,  the  monochromator's  1.25-nm  wide  slits  give  a  50 
A  optical  resolution.  Because  of  its  spectral  scan  from  1.8  to  4.9  pm,  the 
instrument's  resolving  power  R,  defined  as 

R  =  X/AX  ,  (2) 

ranged  from  360  to  980. 

o 

The  50  A  spectral  slit  width  defines  the  width  of  the  theoretical 
triangular  radiant  energy  slit  function  and  is  illustrated  in  Figure  29. 
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Figure  29.  Theoretical  Distribution  of  Radiant  Energ 
with  Wavelength  Exit  Slit 


This  theoretical  triangular  radiant  energy  distribution  indicates  that 

a  peak-normalized  bandwidth  could  be  employed  in  determining  the  flux  measured 

by  the  monochromator.  With  this  peak  normalization  technique  the  effective  50 
o  ... 

A  slit  width  was  centered  about  a  specific  wavelength  which  was  dependent  on 
scan  position.  The  peak  power  which  occurs  at  this  wavelength  was  normalized 
such  that  an  "energy  rectangle's"  height  was  defined.  Hence,  the  instrument's 
spectral  responsivity  R(X)  is  made  equal  to  the  width  of  the  rectangle  R(X) 
enclosing  an  area  equal  to  the  integrated  area  under  the  true  R(x)  area.  This 
normalization  procedure  can  alBO  be  used  when  calculating  the  monochromator's 
calibration  transfer  function. 


The  blackbody's  equation  for  radiance  N..  (T)  is 


r°  2  ci  -2  -1 

na*  <t>  “  -s— -  dX  (Wcm  sr  } 

LX  X5(e  2/XT  -  1) 
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where 


X  i  *  cut-on  wavelength  (cm) 

X  2  a  cut-off  wavelength  (cm) 
ci  =  2irc2h  =  3.741832  X  10“12  (Wcm2) 

c2  =  hc/Kb  -  1.438786  (cmK) 

T  ”  blackboay  source  temperature  (K) 


Equation  3  does  not  have  a  closed  form  solution  and  must  be  solved  through  a 
geometric  series  expansion  before  N^(T)  can  be  input  into  Equation  3  (Ref. 
15).  The  wavelength  interval  between  Xi  and  X2  is  the  monochromator's 
spectral  response  given  earlier.  The  error  introduced  by  bandwidth  normaliza¬ 
tion  is  quantified  when  calculating  the  monochromator's  absolute  responsivity. 


The  normalized  absolute  responsivity  R(X)  over  a  fixed  bandpass  AX  is 

the  ratio  of  the  instrument  output  voltage  to  the  irradiance  E^  or  radiance 

L.,  of  the  calibration  source,  i.  e., 

AX 


R(X)  =  "cal^AXcal 


or 


(4) 


AXcal 

Hence,  the  in-band  irradiance  or  radiance  of  a  source  can  be  computed  from 


EAX source  or  EAX source  *  Vca]./R(X). 
Linearity  is  implicit  in  Equations  4  and  5. 


(5) 


Radiance  from  a  1000  C-rated  Barnes  Model  11-210  blackbody  radiation 
source  with  an  emissivity  e  of  0.99  +  0.01  and  a  0.5-inch  cavity  aperture  was 
used  in  the  absolute  responsivity  calibration. 

Hence,  the  plume  in-band  radiance  relative  to  the  system  transfer 
function  Ks  is  essentially  a  function  representing  the  product  of  the 


normalized  absolute  responsivity  R(a  )  and  its  associated  uncertainties,  shown 
in  Equation  6  with  an  uncertainty  factor  U: 


K8  =  R(X)U  -  I^W/Vd, 


(6) 


where 


N^(T)  =  blackbody  radiance  determined  from  Equation  3;  and 

Vcal  “  calibration  radiance  channel  voltage. 

The  ensuing  system  transfer  function  Ks  is  necessary  for  calibrated  data 
reduction  of  ensuing  plume  radiance  measurements. 

The  uncertainty  factor  U  is  comprised  of: 

(1)  uncertainty  in  determining  the  calibration  constant  Ks,  which 
includes  source  uniformity  and  alignment  uncertainties  in  the  calibration 
system's  optical  transmission  and  inadequate  out-of-band  rejection; 

(2)  effects  due  to  the  unknown  plume  spectral  intensity  curve  (i.  e., 
spectral  peaks) — this  is  largely  mitigated  by  the  linearly  responding 
instrument; 

(3)  radiation  from  outside  the  monochromator's  IFOV  contributing  to 
the  measured  radiance; 

(4)  plume  test  pad  interactions  arising  from  flowfield  impingement 
onto  the  ground  and  reflections  of  radiation  from  the  pad  walls;  and 

(5)  window  contamination  changing  the  calibration  constant. 

The  relative  uncertainties  of  Items  1  and  2  are  discussed  in  detail  in 
Reference  16,  including  a  method  to  calculate  the  effects.  As  mentioned, 
Effect  3  is  minor  with  a  flatter  monochromator  response  curve  R(X).  However, 


the  magnitude  of  Effect  4  is  unknown.  Uncertainty  caused  from  Effect  5  can  be 
reduced  through  frequent  calibration  and  window  cleaning. 


2.4.2  Ultraviolet  (UV)  Instrumentation  -  BATES  plume  UV  spectra  are  of 
interest  both  to  help  verify  solid  propellant  combustion  models,  and  more 
basically,  to  better  construct  a  plume  UV  radiation  model  by  identifying  the 
species  that  radiate  in  the  UV.  Laboratory  UV  emission  spectroscopy  is  being 
undertaken  at  AFRPL  with  selected  solid  rocket  motor  composite  propellants 
(Ref.  17)  and  at  Grumman  Aerospace  Corporation  with  shock-heated  AI2O3 
particles  (Ref.  18).  UV  measurements  have  also  been  taken  with  rocket  motors 
(Refs.  19  and  20).  Reduced  and  minimum  smoke  propellant  motors,  however,  have 
never  been  instrumented  for  UV  measurements  until  this  AFRPL  application. 


Figure  30.  Top  View  of  UV  Measurement  System 

Description  -  The  UV  measurement  system  is  depicted  in  Figure  30  .  It 
is  a  receiving  system  only,  consisting  of  a  Spex  Minimate  Model  1670  grating 
monochromator  and  associated  optics,  and  identical  in  every  way  to  the  IR-E/A 


receiving  unit  except  that  a  BaFl  window  and  a  1200  grooves/mm  diffraction 

o 

grating  blazed  to  operate  with  maximum  efficiency  at  3000  A  were  inserted  in 

place  of  a  fused  silica  window  and  a  300  grooves/mm  grating.  Also,  a  RCA  1P28 

photomultiplier  tube  rated  at  a  spectral  response  exceeding  40m/A/W  from  2600- 
o 

4300  A  was  used  in  lieu  of  the  IR-E/A  systems's  InSb  detector.  Radiation  from 
the  monochromator  exit  slit  was  focused  onto  the  photomultiplier  tube's 
photocathode. 

The  monochromator's  1.25-mm  entrance  and  exit  slits  were  set  for  a  50 
o  o 

A  spectral  resolution.  The  100-ms  2750-4250  A  spectral  scan  was  accomplished 
at  a  rate  of  15  m/s  via  the  same  drive  mechanism  detailed  in  Section  2.4.  The 
difference  in  scan  speed  was  accounted  for  by  a  different  cam  attachment. 

Signal  Conditioning  -  Only  three  voltage  signals  needed  to  be 
recorded:  the  UV  monochromator's  drive  mechanism's  input  voltage,  Vuv^£n;  the 
Vuvx  signal  corresponding  to  the  UV  monochromator  potentiometer's  wiper 
voltage;  and  the  Vuv^  voltage  corresponding  to  the  photomultiplier's  response 
to  the  radiation  focused  on  its  detection  surface.  These  signals  are  sent 
through  300  ft  of  coaxial  cable  to  the  control  room  for  digital  data  sampling 
and  recording  at  250  samples/s  through  a  100  Hz  filter. 

Calibration  -  Absolute  spectral  positioning  was  accomplished  in  a 

manner  similar  to  the  IR-E/A  system  described  in  Section  2.4.3.  VuvUn 

o 

responses  to  3131.55,  3663.28  and  4046.56  A  Hg  line  peaks  were  mapped  vs  Vuvx 
response  outputs,  giving  a  relationship  between  Vuvxin  anc^  Vuv^  Equation 

1,  X  *  a  (Vuvx in/vuvA )  ~  is  used  in  a  linear  curve  fit.  This  established 
an  absolute  spectral  positioning  calibration. 

Figure  31  illustrates  the  dominance  of  first  order  radiation  in  the 
spectral  region  used  here,  eliminating  the  need  for  blocking  filters. 

o 

As  noted  in  the  description  of  the  IR-E/A  system,  the  50  A  spectral 
resolution  defines  the  extent  of  the  theoretical  triangular  radiant  energy 
distribution.  Again,  the  spectral  response  function  R(A )  was  peak  normalized 
and  set  at  a  rectangle  whose  height  was  defined  at  r{A ) .  The  monochromator's 
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Figure  31.  1200  Groovea/mm  Diffracting  Grating 

Efficiency  va  Wavelength 


resolving  power  as  defined  by  Equation  2,  ranged  from  25  and  85  between  2750 
o 

and  4250  A. 


There  is  not  much  radiant  power  from  a  1000  C  blackbody  source  falling 
o 

within  the  2750-4250  A  spectral  range,  but  the  system  calibration  function  Ks 
was  still  calculable  using  the  procedure  outlined  in  Equations  2-6. 

Recognizing  the  potentially  low  signal-to-noise  ratio  with  a  blackbody-based 

* 

absolute  response  calibration,  a  quartz-halogen  calibration  source  was 
originally  proposed.  However,  the  blackbody  source  was  found  to  be 
sufficient. 


The  same  contributions  to  uncertainty  in  determining  the  absolute 
spectral  response  that  were  enumerated  in  Section  2.4.3  also  exist  here. 
Reference  16  describes  how  they  were  incorporated  in  UV  radiance  calibration 
calculations. 


3.0  DATA  REDUCTION  AND  RESULTS 


Data  reduction  methodologies  and  results  are  presented  in  this 
section.  Any  difficulties  in  retrieving  data,  either  from  instrumentation 
problems  or  data  inversion  techniques  are  mentioned  in  this  section  as  they 
strongly  impact  results.  The  four  measurement  areas  presented  here  are: 

(1)  plume  visibility  figure  of  merit; 

(2)  plume  particle  sizing; 

(3)  plume  temperatures  and  partial  pressures;  and 

(4)  plume  UV  signatures. 

3.1  Pluae  Visibility  Figure  of  Merit 

3.1.1  Introduction  and  Data  Reduction  Method  -  One  propellant  ranking 

criterion  is  plume  visibility.  This  is  of  especial  concern  with  low 

visibility  reduced  and  minimum  smoke  solid  rocket  formulations.  In 

conjunction  with  an  ongoing  low  visible  signature  tactical  rocket  motor  grain 
formulation  development  effort,  AFRPL  has  undertaken  solid  rocket  motor 
exhaust  plume  visibility  studies  (Refs.  21  and  22). 

An  experimental  technique  from  which  some  plume  visibility  figures  of 
merit  could  be  deduced  would  greatly  benefit  AFRPL' s  low  visibility  propellant 
development  projects.  The  difficulty  encountered  is  devising  measurements 
which  correlate  well  with  primary  and  secondary  smoke  levels  under  flight 
conditions.  To  address  this  issue  AFRPL  has  developed  an  analytical  computer 
code  aimed  at  evaluating  various  plume  visibility  screening  methods  (Ref. 
23),  Future  efforts  with  radiometric  plume  characterization  measurements  are 
planned. 


It  was  thought  that  the  laser  transmissometer  measurement  could  be 
applied  as  a  plume  visibility  screening  method.  Major  contributors  to  primary 


smoke  visibility  are  the  particle  optical  properties  which  control  the  plume's 
nearfield  contrast.  Contrast  is  a  function  of  scattered  and  absorbed  light. 
The  particle  optical  properties  which  determine  scattering — particle  size,  the 
real  index  of  refraction,  and  the  absorption  coefficient — are  also  the  factors 
which  determine  laser  beam  transmittance.  Moreover,  particle  density  is  an 
affecting  parameter  in  both  plume  visibility  and  laser  transmittance;  these 
two  effects  are  related  and  described  by  Mie  theory  (Appendix  A).  In  fact, 
all  primary  smoke  plume  visibility  factors  except  afterburning  mechanisms 
relate  to  laser  beam  transmittance. 

There  are  major  drawbacks  in  applying  the  laser  transmissometer 
experiment  in  this  way.  Secondary  smoke  arising  from  water  and  acid 
condensation  in  the  plume  farfield  is  not  addressed,  nor  is  the  presence  of 
ignition  smoke  or  plume  scattering  from  either  direct  or  diffuse  sunlight 
considered.  Also,  the  firings  under  this  project  are  static;  flight 
conditions,  which  have  a  major  impact  on  the  plume  flowfield  and  therefore  the 
visible  signature,  are  not  accounted  for.  The  laser  transmission  is  made  only 
from  the  90°  aspect  angle,  hence  ignoring  other  radiative  flux  directions 
which  contribute  to  total  plume  visibility  (such  as  the  rocket  chamber 
"searchlight"  effect  and  the  very  real  angular  dependence  of  light 
scattering).  Additionally,  laser  transmission  is  monochromatic,  and  although 
the  particle  optical  properties  on  which  transmittance  is  dependent  apply  to 
other  visible  wavelengths,  the  mathematical  scattering  efficiency  functions 
describing  Mie  theory  are  transcendental  and  highly  sensitive  to  wavelength 
as  well  as  particle  size  and  aspect  angle.  A  final  point  which  is  actually  a 
combination  of  the  drawbacks  listed  above  is  that  any  correlation  between 
laser  beam  extinction,  its  application  as  a  secondary  measure  of  primary  smoke 
to  the  ultimate  plume  visibility  criterion,  total  plume  detectability,  has  not 
been  undertaken.  This  is  the  bottom  line  which  a  plume  visibility  screening 
figure  of  merit  must  address. 

It  was  decided  to  apply  the  laser  transraittances  from  the  AFRPL 
particle  size  measurements  to  assess  plume  visibility  with  the  above  caveats 
taken  into  account.  Laser  transmissometer  instruments  and  geometry  were 
detailed  in  Section  2.3.1  of  this  report. 
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3.1.2  Transmissometer  Results  and  Interpretation  -  Laser  transmission  data 
was  successfully  collected  for  29  firings.  In  each  of  these  firings  the  laser 
action  period  was  divided  into  ten  intervals,  excluding  the  time  interval  the 
shutter  breaker  blocked  the  laser  output.  The  mean  transmissivities  of  these 
intervals  and  the  one  sigma  values  are  recorded  in  Table  5  as  well  as  the  mean 
transmissivities  and  standard  deviations  within  each  propellant  group.  Note 
that  while  these  mean  transmissivity  numbers  may  reflect  accurate  figures  of 
merit  for  a  large  ensemble,  the  small  samples  and  the  variations  of  nozzle 
geometries  here  detract  from  a  straightforward  comparison  of  different  types 
of  propellants.  The  uncertainties  are  reflected  in  the  propellant  series' 
standard  deviations  but  only  to  the  extent  of  sample  size  and  only  in  cases 
where  sample  size  exceeds  one.  It  may  not  be  possible  to  ascribe  the  one 
sigma  statistical  description  of  the  mean  transmissivities  to  a  visibility 
criterion  because  all  sample  sizes  are  small. 

Considering  the  just  mentioned  caveat,  interesting  transmissivity  (and 
therefore,  plume  visibility)  comparisons  can  still  be  made.  Table  5  lists 
three  reduced  smoke  sidewinder  propellants.  Reviewing  each  propellant 
formulation's  mean  transmissivity  xp  and  the  data  quality  privisos,  both  AFRPL 
RSS  (tp  =  83.054)  propellant  plumes  are  attenuated  to  a  greater  extent  with 
our  Ar+  laser  transmissometer  0.5145  pm  than  is  Aerojet's  RSSA  formulation. 
By  this  measure  only,  it  appears  Aerojet's  RSSA  formulation's  plume  is  less 
visible. 

Similarly,  of  the  propellant  formulations  tested  without  solid 
additives,  the  MSMP-4  (xp  =  96.474)  and  DWP  (xp  °  96.429)  formulations  proved 
superior  in  our  visibility  figure  of  merit  than  did  KKN  (xp  ■  92.389)  or  RS 
GAP  80  ( xp  =  90.538).  It  must  be  noted  that  the  RS  GAP  80  transmissometer 
experiment  was  configured  with  the  laser  beam  path  2  cm  (instead  of  2  m)  from 
and  parallel  to  the  nozzle  exit  plane. 

It  is  important  to  once  again  point  out  that  a  true  statistical 
description  of  representative  plume  laser  transmission  for  each  of  the 
propellant*-  types  tested  was  not  accomplished  with  the  limited  data  presented 
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Table  5.  Tranaai.88iviti.es  for  29  15-Pound  BATES  Motor  Firings 


AFRPL  PROPELLANT  MEAN 

MOTOR  TRANSMISSIVITY 


IDENTIFICATION 

PROPELLANT 

TRANSMISSIVITY  (Zx) 

°T 

Z  (Xjj) 

°T 

144B-008 

*RSS 

83.049 

0.88 

83.049 

179B-001 

*RSSA 

88.970 

1.08 

91.695 

3.68 

179B-004 

*RSSA 

94.169 

1.03 

91.695 

3.68 

194B-001 

*ARCOCEL-423 

87.157 

5.40 

84.054 

4.39 

194B-002 

♦ARCOCEL-423 

80.951 

1.74 

84.054 

4.39 

240B-004 

CWB  (DTS-9886) 

85.673 

2.34 

85.673 

— 

241B-001 

KKN  (DTS-9889) 

94.968 

0.750 

92.389 

2.89 

241B-002 

KKN  (DTS-9889) 

95.886 

1.35 

92.389 

2.89 

241B-004 

KKN  (DTS-9889) 

89.371 

2.04 

92.389 

2.89 

241B-005 

KKN  (DTS-9889) 

90.281 

1.32 

92.389 

2.89 

241B-006 

KKN  ( DTS-9889 ( 

91.441 

0.830 

92.389 

2.89 

242B-001 

DWP  (TPH-8289) 

95.573 

1.76 

96.429 

1.73 

242B-002 

DWP  (TPH-8289) 

94.710 

1.27 

96.429 

1.73 

242B-003 

DWP  (TPH-8289) 

99.194 

0.700 

96.429 

1.73 

242B-005 

DWP  (TPH-8289) 

93.039 

3.00 

96.429 

1.73 

242B-007 

DWP  (TPH-8289) 

96.571 

1.24 

96.429 

1.73 

242B-008 

DWP  (TPH-8289) 

99.179 

0.843 

96.429 

1.73 

242B-009 

DWP  (TPH-8289) 

95.728 

1.99 

26.429 

1.73 

242B-010 

DWP  (TPH-8289) 

97.514 

1.78 

96.429 

1.73 

242B-011 

DWP  (TPH-8289) 

96.389 

2.68 

96.429 

1.73 

242B-013 

DWP  (TPH-8289) 

95.315 

3.95 

96.429 

1.73 

242B-016 

DWP  (TPH-8289) 

94.334 

1.01 

96.429 

1.73 

242B-018 

DWP  (TPH-8289) 

98.158 

0.432 

96.429 

1.73 

242B-019 

DWP  (TPH-8289) 

95.123 

2.36 

96.429 

1.73 

242B-020 

DWP  (TPH-8289) 

95.123 

2.36 

96.429 

1.73 

242B-022 

DWP  (TPH-8389) 

97.728 

1.09 

96.429 

1.73 

242B-023 

DWP  (TPH-8289) 

97.562 

1.26 

96.429 

1.73 

248B-001 

MSMP-4 

96.474 

0.704 

96.474 

— 

262B-005 

RS  GAP  80 

90.538 

3.39 

90.538 

— 

♦These  are  reduced  smoke  sidewinder  propellants.  Firing  number  144B-008  was 
an  AFRPL  formulation  with  0.5  weight  percent  A1  additive  loading;  the  179B- 
series  was  an  Aerojet  formulation  with  0.5  weight  percent  ZrC;  and  the  194B- 
series  was  an  ARC  formulation  with  1.0  weight  percent  ZrC  (Ref.  1). 


here;  hence,  the  extent  of  error  in  the  above  observations  cannot  yet  be 
reported.  In  addition,  the  link  between  plume  transmissivity  and  an  accurate 
visibility  figure  of  merit  option  is  not  established. 


One  final  point:  previous  efforts  to  accurately  correlate  the  mean 

particle  diameter  to  such  motor  and  performance  parameters  as  nozzle  throat 
diameter,  particulate  concentration  within  the  motor  chamber,  average  chamber 
pressure  and  average  chamber  residence  time  (Ref.  24)  arose  from  observing 
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that:  a  reasonable  single-variate  correlation  between  the  particle  and  nozzle 
throat  diameters  may  exist  (Ref.  25).  This  leads  to  the  possibility  that  a 
relationship  between  transmissivity  and  these  parameters  may  exist.  Data 
scatter,  however,  was  the  limiting  factor.  On  second  glance  some  data  scatter 
would  appear  reasonable.  The  relationship  between  laser  transmissivity  and 
mean  particle  diameter  is  not  only  logarithmic  but  also  depends  on  the 
difficult  to  solve  Mie  theory  equations.  This  will  be  shown  in  Section  3.2.1. 

3.2  Plume  Particle  Sizes 

3.2.1  Transmissometer  Data  Reduction  Methods  -  The  laser  transmission/ 
scattering  measurement  described  in  Section  2.3.1  present  the  in-situ  particle 
sizing  data  described  in  this  section.  Two  data  reduction  schemes  were 
envisioned  to  deconvolve  the  measured  angular  and  transmitted  light 
intensities  to  particle  sizes:  an  AEDC-supplied  computer  code  designed  to 

yield  PSDF  histograms  and  an  Aerospace  Corporation-contracted  computer  code 
designed  to  yield  mean  particle  sizes.  For  reasons  to  be  shown  in  this 
section  only  the  second  of  these  techniques,  the  Aerospace  Corporation  method, 
was  used  to  retrieve  particle  sizes. 

The  AEDC-supplied  size  deconvolution  code  (SIZFRED)  (Ref.  26)  input 
parameter  requirements  were  used  to  help  design  the  AFRPL  laser  scattering 
configuration.  SIZFRED  was  generically  designed  to  yield  PSDF  bins  from 
angularly  scattered  and  orthogonally  polarized  light  incident  on  a  multiangled 
array  of  signal  detection  channels,  all  normalized  by  the  radiation  intensity 
from  an  unpolarized  90°  scattering  channel.  The  laser  light  was  to  be 
initially  oriented  45°  to  the  scattering  plane  so  both  parallel  and 
perpendicular  components  would  be  present,  thus  simulating  unpolarized  light. 
The  scattering  volume  wa9  to  be  narrowly  defined  by  a  focused  array  of 
detectors.  The  AFRPL  scattering  measurement  apparatus  was  modeled  in  this 
way. 

SIZFRED  was  later  modified  to  two  families  of  computer  codes: 
(l)  MIESTEP,  which  calculates  Mie  functions  for  the  polarization  states  given 
the  detection  angular  positions  and  FOV's;  and  (2)  SIZIM,  which  performs  the 
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PSDF  deconvolution  from  the  MIESTEP  data  tape  through  a  modified  Phillips- 
Twomey  (Ref.  27)  constrained  linear  inversion  technique. 

Although  the  two  codes  are  described  in  some  detail  by  Curry  (Ref. 
27),  explanations  of  Mie  scattering  theory  and  the  constrained  linear  in¬ 
version  numerical  method  are  in  order  to  explain  the  applicability  of  the  AEDC 
deconvolution  methods  to  the  AFRPL  laser  scattering  experiment.  The  tenants 
of  Mie  scattering  as  they  specifically  apply  to  the  AFRPL  plume  particle 
sizing  effort  are  given  in  Appendix  A. 

To  understand  the  constrained  linear  inversion  technique  one  must 
recognize  that  the  scattered  intensity  components  can  be  written  as  Fredholm 
integral  equations  in  the  form  (Ref.  28) 

Gi(yi)  =  JQK(yi,  D)dD.  (7) 

In  this  specific  particle  sizing  application,  the  scattering  intensity 
components  for  each  scattering  detector  are  Fredholm  integral  equations  of 
like  form: 

IX(e)  *  c  jIaPj_(0)n(a)da  (8) 

and 

I„(6)  ■  c[oapHn^a^  da  (9) 

where 

6  =  scattering  detector  aspect  angle, 

C  »  scattering  detector  absolute  spectral  response  calibration 

constant, 

Op  ,  (6 )  ■  differential  scattering  cross  sections; 


a 


dimensionless  size  parameter 


a  -  2irr/X  where  r  is  particle  radius;  and 
n(a)  =  particle  size  distribution  function. 


An  explanation  of  the  terms  6,  op,  a  and  n(a)  is  contained  in  Appendix  A. 


The  constrained  linear  inversion  technique  embodied  in  SIZIM 
essentially  begins  with  the  polarization  intensity  ratio  definition: 


P(6) 


IA(0)  -  I,i  (0) 

FJTej  +  i  ,|  (e) 


do) 


Substituting  Equations  8  and  9  into  Equation  10  yields  the  relation, 

C  aPl(9)n(a)d  -  f“  ap  (0)n(a)da 

P(0)=7_b - i - — -  (H) 

Jo  aPj.(0)n{a)da  “Jo  Op;|(6)n  (“Jda 

which  explicitly  relates  measured  data  P(0)  to  the  PSDF  n(a).  The  PSDP  is 
retrieved  from  the  equation  through  the  Phillips-Twomey  constrained  linear 
inversion  method.  The  method  is  modified  to  provide  a  linear  solution  for  the 
nonlinear  Equation  11.  The  PSDF  n(a)  is  retrieved  by  first  holding  the 
numerator  constant  and  solving  for  n(o ) ,  then  holding  the  denominator  constant 
using  the  just-retrieved  n(a )  as  an  iteration  starting  point.  This  recursive 
approach  is  terminated  upon  convergence  to  a  solution. 


Young  (Ref.  29)  has  approached  what  is  essentially  the  solution  of 
Equation  11  differently.  He  has  approximated  the  nonlinear  integral  equation 
by  using  a  first-order  linear  expansion  and  has  shown  that  the  Fredholm  kernal 
(given  in  Equation  11  as  op^(e)  and  Op  ^(e))  varies  more  than  Curry  (Ref.  26) 
has  determined.  Young  attributes  this  difference  to  his  not  making  the 
assumption  of  holding  the  denominator,  then  numerator,  constant  in  successive 
recursions  as  Curry  does.  The  impact  of  this  analysis  on  the  AFRPL  scattering 
experiment  is  to  shed  uncertainty  on  the  original  particle  size  retrieval 
method  of  choice — the  MIESTEP/SIZIM  family  of  codes. 


In  his  analysis,  Young  ran  uncalibrated  relative  intensity  synthetic 
data  like  that  acquired  in  the  AFRPL  laser  scattering  set-up  through  his 
first-order  linearized  particle  size  inversion  scheme  to  determine  if  a 
relative  size  distribution  could  be  accurately  retrieved.  His  results,  shown 
in  Figure  32,  point  to  the  fact  that  in  this  case,  data  from  the  laser 
scattering  experiment  yields  size  distribution  information  only  for  particles 
in  the  0.3  to  0.5  m  size  range.  This  is  very  little  information  indeed. 
Moreover,  the  data  used  in  this  analysis  was  noiseless,  so  error  propagation 
through  the  linear  deconvolution  algorithm  was  not  taken  into  account. 


The  perfect  noiseless  data  necessary  for  the  aforementioned  inversion 
was  not  supplied  by  the  AFRPL  experiment.  Figure  33  depicts  the  low  signal- 
to-noise  ratio,  showing  very  little  difference  in  scattering  detector  pre-  and 
post-ignition  signal  levels.  Appendix  C  addresses  the  AFRPL  efforts  in 
enhancing  scattering  detector  signals  through  an  innovative  adaptive  noise 
cancellation  scheme. 
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Figure  33 .  Scattering  Detector  vs  Time  Signal  with  a  15-lb 
BATES  Motor  Chamber  Pressure  Trace 


Given  the  debilitating  uncertainty  within  the  PSDF  inversion  scheme 
combined  with  the  experimental  uncertainties  mainly  from  the  low  scattering 
detector  signal-to-noise  ratios  (but  also  from  unknowns  in  the  calibration 
procedure  and  scattering  focal  volume  determination),  it  was  decided  to  use 
the  laser  scattering/transmissometer  experiment  to  yield  a  volume-to-surface 
(or  Sauter)  mean  particle  diameter  integrated  over  the  laser  line  of  sight 
through  the  plume  (Refs.  29  and  30). 


The  single-color  transmissometer  algorithm  described  in  these  refer¬ 
ences  has  been  codified  into  a  mean  particle  size  retrieval  computer  program, 
A3 2 CODE  (Ref.  31). 


The  basic  retrieval  method  used  in  A32CODE  recognizes  the  low  mass 
loading  conditions  prevalent  in  the  low-visibility  propellant  two-phase  plume 
flows.  These  conditions  make  reasonable  the  application  of  single-scattering 
from  homogeneous,  spherical  particle  assumptions  (discussed  in  Appendix  A)  and 
the  application  of  Beer's  Law  toward  describing  laser  light  extinction: 


lrrr 


Qext  (X»r)  n(r)dr 


(12) 


where 


t  ■  transmittance; 


L  ■  laser  path  length; 


r  *  particle  radius; 

Qext  *  extinction  efficiency; 
A  *  laser  wavelength;  and 
n(r)  *  PSDF. 


Salient  concepts  behind  Qexj  And  n(r)  are  discussed  in  Appendix  A. 


In  Young' 8  methodology  (Ref.  29),  he  defines  the  mean  extinction 
efficiency  Qext^ )  as  being  a  weighted  mean  value: 

- .  .  Jo'2  5e«,A-  cl"(rHr 


Qext(*) 


r  n(r)dr 


The  total  mass  loading  is  the  intdgratdd  product  of  thd  characteristic  bulk 
density  of  the  particulate  material  (Pp),  the  volume  of  the  particle  and  the 
number  of  particles  is  expressed  as 


52 


Hence,  A32CODE  has  as  its  input  parameters  (Ref.  32):  (1)  plume 
diameter  L  at  laser  transmissometer  station  (cm);  (2)  particulate  mass  loading 
(g/crn-*);  (3)  particulate  material  bulk  density,  pp  (g/cm^);  and  (4)  laser 
transmittance  through  the  plume.  Additional  inputs  include  real  and  imaginary 
components  of  the  particulate  index  of  refraction,  laser  wavelength  and  choice 
of  unimodal  or  bimodal  rectangular  size  distributions.  Uncertainties  in 
measurement  are  imbedded  in  error  base  input  options  for  Parameters  1-4,  and 
lack  of  knowledge  of  indices  of  refraction  can  also  be  specified  in  the  code 
inputs.  The  final  output,  mean  particle  radius  I32’  reflects  the  input  error 
budgets. 


It  should  be  cautioned  again  that  the  Beer's  law  assumptions  in  the 
physics  behind  the  code — a  low  concentration  plume,  a  homogeneous  plume  media 
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(with  an  index  of  refraction  n“l-0k)  and  the  interaction  mechanisms  consisting 
of  absorption  and  scattering — limit  the  applicability  of  the  code  to  the  low- 
visibility  type  propellants  used  in  this  study.  Furthermore,  the  single 
scattering  by  homogeneous,  solid  spherical  particles  implied  by  using  a  mean 
extinction  efficiency  Qext^ )  introduces  enough  uncertainty  that  results  must 
be  couched  in  terms  of  effective  mean  volume-to-surface  particle  radius  over 
the  laser  line  of  sight.  These  objections  aside,  the  just  described  analytic 
procedure  is  sound  and  helps  provide  AFRPL  with  its  first  substantive  in-situ 
particle  sizing  capability,  from  measurement  to  final  data  reduction. 

3.2.2  Transmissometer  Results  and  Interpretation  -  Considerable  difficulty 
was  experienced  in  taking  suitable  BATES  motor  laser  transmission  data.  With 
solid  Al  or  ZrC  additive  loadings  approaching  5%  by  weight,  laser 
transmissivity  was  virtually  zero,  rendering  this  type  of  data  unusable. 
Since  a  significant  number  of  AFRPL  BATES  firings  are  not  the  solid  loaded  low 
smoke  type,  many  targets  of  opportunity  were  necessarily  bypassed.  Also,  the 
laser  transmissometer  configuration  described  in  Section  2.3  worked  only  for 
15-lb  BATES  motors.  Actions  taken  in  1984  to  make  the  transmissometer  axially 
adjustable  for  30-  and  70-lb  BATES  motors  were  effected  late  enough  to  not 
affect  any  solid  additive  loaded  plume  measurements  reported  here.  As  a 
consequence,  of  the  29  firings  listed  earlier  in  Table  5  (Section  3.1.2)  with 
usable  laser  transmission  data,  only  five  had  Al  or  ZrC  additives;  the  other 
24  had  no  solid  propellant  ingredients.  Table  6  summarizes  the  particulate 
sizes  derived  from  the  A32CODE  described  in  the  previous  section. 

In  reducing  the  data,  input  error  bounds  were  estimated. 
Transmissivity  standard  deviations  served  as  the  error  bands  for  those 
measurements.  Additionally,  particulate  mass  loadings  were  determined  by 
elemental  molecular  weights  and  gas  and  particulate  matter  exhaust  densities 
were  calculated  from  an  AFRPL  theoretical  one-dimensional  isentropic 
equilibrium  program  (Ref.  32).  A  llOZ  error  range  was  assigned  to  this 
computation  and  to  specific  gravities  for  pure  form  AI2O3  and  Zr02»  Indices 
of  refraction  ranging  between  (1.7-1.8)n  +  i  (0-0.02)k  were  assigned  to  the 
particles.  Finally,  the  length  of  the  laser  path  through  the  plume  was 
determined  to  be  9.75  1.  0.65  cm. 
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Table  6.  Particulate  Size  fron  Transmissometer  Ezperiaent 


PROPELLANT 

EXHAUST 

COMPUTED 

AFRPL  MOTOR 

PARTICLE 

PARTICULATE 

D32  (ym)  OF 

IDENTIFICATION 

LOADING 

PRODUCT 

EXHAUST  PARTICLES 

144B-008 

0.502  A1 

AI2O3 

0.700  +  0. 106(+15 .12) 

179B-001 

0.502  ZrC 

Zr02 

0.744  +  0.139(^18.82) 

179B-004 

0.502  ZrC 

Zr02 

2.08  +  0.56  (+  27.12) 

194B-001 

1.02  ZrC 

Zr02 

1.63  +  ,084(+  51.72) 

194B-002 

1.02  ZrC 

Zr02 

0.922  +  0.134(+29.02) 

It  must  be  remembered  that  these  particle  size  measurements  were  taken 
from  a  thus  far  unverified  laser  transmissometer  technique.  Implicit  in  this 
caveat  are  the  assumptions,  both  analytical  and  physical,  described  in  Section 
3.1  and  Appendix  A.  Although  we  are  confident  of  the  size  retrieval 
methodology,  another  corresponding  measurement  or  absolute  system  calibration 
has  yet  to  be  adequately  done. 

The  small  sample  size  of  five  coupled  with  three  different  values  for 
the  propellant  particle  loading  variable  make  any  trend  analysis  of  the 
derived  particle  size  difficult  and  highly  speculative 

Table  7  lists  the  comparison  between  the  transmissometer  measurement- 
derived  mean  particle  diameter  and  the  semi-empirical  correlation  diameter 
calculated  from  Equation  16.  The  above  listed  mass-volume  diameter  D43  was 
chosen  by  Nickerson  (Ref.  24)  instead  of  the  volume-surface  diameter  D32 
because  it  serves  to  better  predict  two-phase  flow  losses,  but  the  difference 
should  not  be  as  pronounced  as  shown  above.  The  D32  values  exceed  the  D43 
diameters  by  several  factors  in  each  of  the  three  firings.  The  question  of 
which  figure  is  better  representative  of  the  true  rocket  exhaust  particle  size 
cannot  yet  be  inarguably  answered.  Given  this  limited  information  no  attempt 
was  made  to  correlate  the  computer  mean  particle  size  with  motor  geometry  or 
performance  parameters.  Nickerson  derived  such  an  empirical  relationship  used 
in  the  above  tabulation  with  66  particulate  samples  captured  in  the  farfield 
from  a  variety  of  solid  rocket  motor  exhaust  flows;  this  curve-fit  is  embedded 
in  the  two-dimensional /two-phase  Solid  Performance  Program  (SPP)  used  to 
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Table  7.  D32  ( from  Transmittance)  and  D43  (from  Nickerson)  Comparison 


AFRPL  MOTOR 

IDENTIFICATION 

EXHAUST  PARTICULATE 

PRODUCT 

D32  (pm) 

(TRANSMISSOMETER) 

D43  (pm) 
(CORRELATION) 

144B-008 

AI2O3 

0.700 

0.293 

179B-001 

Zr02 

0.744 

0.222 

179B-004 

Zr02 

2.08 

0.0988 

194B-001 

Zr02 

1.63 

0.280 

194B-002 

Zr02 

0.922 

0.276 

to  predict  motor 

performance  and  an  input 

to  JANNAF  flowfield 

and  radiation 

prediction  computer  codes.  Nickerson's  relation  is: 

D43  -  3.6304  »t0-2532  <1  -  .-0-00 <»»««  c  PC,  )  (16) 

where 

D43  **  mass-weighted  average  diameter  (pm); 

Dt  =*  nozzle  throat  diameter  (in); 

5C  *  AI2O3  or  Zr(>2  concentration  in  chamber  (g-mol/100  g); 

Pc  =*  average  chamber  pressure  (psia);  and 

t  *  average  chamber  residence  time  (ms). 

Particle  capturing  techniques  used  in  establishing  data  points  to  derive 
Equation  16  suffer  from  the  same  deficiencies  enumerated  in  Section  2.3  (and 
in  fact,  may  suffer  from  even  more  bidses  because  some  particles  were  gathered 
from  petri  dishes  lying  exposed  in  the  plume  flow).  The  accurate  specific 
impulse  predictions  of  SPP  point  to  the  semi-empirical  model's  authenticity 
but  only  toward  measuring  specific  impulse.  More  data  points  must  be  acquired 
to  truly  verify  Equation  16.  Even  this,  however,  is  no  substitute  for  the 
necessary  step  of  taking  an  accurate  calibration  or  parallel  measurement  to 


gain  confidence  in  the  laser  transmissometer  experiment's  mean  particle 
diameter  results. 


3.2.3  Particle  Collection  Probe  Technique  Data  Reduction  Methods  -  The 
subsonic  particle  collection  probe  described  in  Section  2.3.2  captured 
particulate  samples  for  AFRPL  in-house  analysis.  At  the  time  of  this  writing, 
the  AFRPL  image  analysis  capability  could  not  give  a  size  distribution  count 
over  the  full  extent  of  the  expected  size  range.  This  is  unfortunate  because 
correlating  mean  sizes  with  the  A32C0DE  effective  mean  size  determination  of 
the  previous  section  would  enhance  confidence  in  the  laser  transmissometer 
measuring  technique's  results.  The  data  would  be  useful  despite  biases 
introduced  by  dust  entrainment,  agglomeration  and  bow  shock  effects.  To 
overcome  these  error  sources,  a  sharp-edged  supersonic  particle  prove  as 
described  in  Section  2.3.2  is  currently  in  fabrication  and  will  be  applied  at 
AFRPL  in  the  future  along  with  an  upgraded  image  analyzer. 

The  possible  deposition  of  carbon,  aluminum  chloride  or  other  surface 
contaminants  on  the  particles  was  not  analyzed  via  any  x-ray  electron  probe 
technique,  as  AFRPL  does  not  presently  possess  the  capability. 


3.2.4  Particle  Collection  Probe  Technique  Results  and  Interpretation  - 


Photomicrographs  of  captured  rocket  exhaust  particles  indicate  that  a  few  may 
be  hollow  (Fig.  33).  Dawbarn  and  Kinslow  (Ref.  33)  noted  the  same  phenomenum 
in  their  work  and  have  summarized  several  suspected  mechanisms.  Alumina  has 
been  postulated  to  form  a  stable  molten  core  of  aluminum  on  the  surface, 
expanding  by  pressure  from  aluminum  vapor;  this  appears  to  be  unlikely  due  to 
the  severe  aerodynamic  forces  interacting  with  particulates  within  the  exhaust 
plume  flowfield.  A  more  likely  explanation  is  that  the  AI2O3  inflation  occurs 
from  gases,  H2  or  CO2,  within  the  plume.  These  inflated  AI2O3  particles  could 
be  nonstoichiometric  a  -phase  alumina  not  reacting  with  these  gases. 
Nonetheless,  inherent  in  all  these  suggestions  is  that  hollow  spheres  arise 
from  aluminum  burning  outside  the  combustion  chamber,  therefore  accounting  for 
but  a  fraction  of  the  total  amount  of  alumina  captured.  This  seems 
reasonable — hot  particle  streaks  in  solid  rocket  motor  plumes  is  well- 
observed  phenomena  at  AFRPL  and  elsewhere. 


.  well- 


Figure  34.  SEM  Photomicrograph  of  Two  Captured  10-um  Al^O^  Spheres. 

Note  that  the  one  on  the  right  appears  to  be  hollow. 

* 

Reference  34  cites  high  ambient  humidity  (greater  than  80%)  as  a 

strongly  contributory  factor  in  the  formation  of  hollow  AI2O3  spheres,  since 
they  observe  that  aluminum  spheres  burn  better — and  hence  more  hollow  spheres 
are  produced — in  a  wet  atmosphere.  The  AFRPL  findings  dispute  this.  All  data 
reported  here  were  taken  at  relative  humidity  levels  usually  well  below  30%. 
Also  the  high  H2O  content  of  a  rocket  plume  logically  disputes  the  empirical 
high  relative  humidity-hollow  particle  correlation;  the  effective  and  already 
high  humidity  within  the  plume  should  be  relatively  independent  of  the 

external  environment. 

Ramifications  of  hollow  particles  with  a  cavity  opening  such  as  shown 
in  Figure  34  are  substantial.  The  effective  emissivity  of  an  aggregation  of 
AI2O3  particles  with  this  characteristic  would  be  much  higher  than  pure  form 

AI2O3,  as  the  cavity  could  emit  blackbody  radiation.  The  hole  also  makes  the 

particles  obviously  nonspher ica 1 ,  and  Mie  theory  mechanisms 


described  in  Appendix  A  might  not  adequately  explain  plume  scattering. 
Electromagnetic  radiation  may  be  absorbed  into  the  cavity,  then  reradiated 
thermally;  also,  it  may  be  reflected  several  times  within  the  cavity,  exiting 
at  different  angles  than  postulated  by  scattering  theory.  Additionally,  mass 
conservation  assumptions  in  flowfield  and  radiation  prediction  computations 
(i.  e.,  differences  in  predicted  particle  sizes  between  hollow  and  solid 
spheres)  makes  the  presence  of  hollow  spheres  a  significant  phenomenon. 
Assumptions  regarding  the  particle  density  impact  not  only  on  results  from 
radiation  predictions  but  also  on  the  A32  code  reported  earlier  in  this 
section.  Strand  (Ref.  34)  reports  densities  of  hollow  AI2O3  rocket  exhaust 
particles  which  vary  from  1.5  to  3.5  g/cm^,  differing  substantially  from  the 
3.7  g/cm^  value  used  in  A32  code  to  determine  the  mean  particulate  size. 

3.3  Plume  Temperatures  and  Particle  Pressures 

3.3.1  IR-E/A  Technique  Data  Reduction  Methods  -  As  described  in  Section  2.4 
of  this  report,  an  IR-E/A  fixed  line-of-sight  multispectral  measurement 
technique  was  used  during  these  firings.  The  spectral  emission/absorption 
scan  data  was  to  be  deconvolved  to  temperature  and  partial  pressure  profiles 
through  an  iterative  Abel  inversion  algorithm  developed  by  Young  (Ref.  35). 

Young's  methodology  can  be  illustrated  by  the  gas-only  radiance  LP  and 
transmittance  xP  relations  developed  in  Appendix  B.  These  relations  can  be 


dxP( s) 


expressed  as: 


L  A(A’  T  }  =  J  (LBBa(s>  +  l\(s))  ds 


where 


^  =  wavelength; 


Tp  =  plume  temperature; 


s  =  line  of  sight; 


1  =  distance  along  s; 


Also,  in  tP<s)  -  -f  cp(s)  pp(s)  Kp( s)  yp(s)d(s),  (18) 

^  o 

where 

cP(s)  =  concentration  at  a; 
pp(s)  “  pressure  at  s; 

Kp(s)  *  monochromatic  absorption  coefficient  at  s;  and 

yP(  s)  =  monochromatic  band  model  function  at  s. 

Young's  algorithm  solves  Equation  18  for  pp(s)  given  a  known  concentration 
profile  cp(s)  and  can  henceforth  also  be  used,  as  indicated  in  Appendix  B,  to 
retrieve  temperature  profiles.  The  Abel  inversion  scheme  described  in 
Reference  35  essentially  uses  the  -d  xP(s)/ds  term  in  Equation  16  as  the 
weighting  function  in  the  inversion  procedure.  Note  that  xp  is  written  as  a 
function  of  both  wavelength  X  and  position  s  along  the  plume  line  of  sight. 

Young  (Ref.  35)  calculates  the  near  infrared  bandpasses  at  which  the 
major  plume  constituents  H2O,  CO2  and  CO  radiate  without  overlap  at  the 
elevated  temperatures  and  pressures  expected  with  BATES  motors  as: 

(1)  H20:  2.3-2.64  pm 

3. 1-3.5  ym 
6. 0-6. 9  ym 

(2)  C02:  4.14-4.23  ym 
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(3)  CO: 


5. 0-5. 3  ym. 


These  define  the  spectral  regions  under  which  an  IR-E/A  inversion  can  be 
accomplished,  with  single  plume  species  radiators  and  therefore  enough  of  a 
pure  -d  t^(s)/d8  weighting  function  variation  along  X  to  justify  a 
multispectral  experimental  technique. 

Young  also  models  a  radial  temperature  profile  along  s  with  the  plume 
centerline  peak  temperature  at  1500K  and  declining  monotonically  to  750K  at  a 
10-cm  radius.  He  also  provides  synthetic  IR-E/A  radiance  data  for  6.3  ym  H2O, 
4.3  ym  CO2.  Despite  this  reasonable  BATES  temperature  profile  and  radiance 
model,  the  ensuing  weighting  function  -di^sj/ds  does  not  vary  enough  to 
extract  accurate  temperature  or  concentration  profiles,  as  illustrated  in 
Figures  35  and  36.  There  is  a  1302  error  bound  even  in  this  best  case.  The 
less  stressing  nonscanning  monochromatic  case  was  then  analyzed,  and  a  similar 


Figure  35.  Temperature  Inversion  Results  from  Young's 
BATES  Plume  Model  (Ref.  25) 
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lack  of  -dxP(s)/ds  resolution  leads  to  only  slightly  better  inversion 
accuracies.  These  poor  inversion  results  occurred  even  with  noise-free 
synthetic  data.  The  plume  analyzed  here  was  particle-free  and  represents  a 
simpler  model  that  one  with  particle  emission  and  absorption  effects. 

The  conclusion,  which  was  reached  midway  through  the  AFRPL 
experimental  efforts,  was  that  a  multispectral  IR-E/A  temperature  and  partial 
pressure  plume  profile  could  not  be  accurately  retrieved  with  solid  rocket 
motors.  The  alternatives  to  pursue  were  essentially  two:  convert  the  IR/E-A 
apparatus  to  multiple  in-band  and  spatially  scanning  emission  and  absorption 
measuring  detectors  similar  to  that  applied  at  AEDC  (Ref.  36)  or  abandon  the 
AFRPL  IR-E/A  efforts. 


Figure  36.  Concentration  Inversion  Results  from  Young's 
BATES  Plume  Model  (Ref.  26). 


AFRPL  chose  to  pursue  the  analytic  portion  of  a  gas/particle 
temperature  and  concentration  profile  retrieval  method  despite  probable  new 
instrumentation  requirements.  The  analysis  (Ref.  37)  and  computer  code  (Ref. 
38)  centered  around  a  spatially  scanning  IR-E/A  measurement  scheme.  The  AFRFL 
experimental  apparatus,  however,  remained  the  same  due  to  cost  and  space 
constraints.  The  AFRPL  scanning  monochromator  was  locked  at  the  2.7-pm  p 
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branch  H2O  band  and  used  as  a  narrow  bandpass  radiometer  to  gauge  the  water- 
band  intensity  variation  through  the  BATES  firing. 


3.3.2  IR-E/A  Technique  Results  and  Interpretation  -  Because  the  IR-E/A 
technique  could  not  be  used  as  originally  expected  to  yield  plume  temperatures 
and  partial  pressures,  no  data  of  this  type  can  be  presented.  However,  some 
infrared  emission  data  in  the  2.7-pm  H2O  band  was  taken  by  the  blackbody 
source/chopper/ lock-in  amplifier  arrangement  described  in  Section  2.4.1. 

There  was  difficulty  in  assigning  system  calibration  factors  to  the 
monochromator's  voltage  output  because  of  inconsistent  wavelength  and  response 
calibration  techniques.  Hence,  a  quantitative  firing-by-firing  comparison  is 
not  possible,  especially  when  considering  that  the  plume  station  radiation 
depends  on  variable  motor  ballistic  and  geometry  parameters.  The  predictions 
of  Figure  37  illustrate  how  a  2.4-2. 7  pm  bandpass  station  radiation  can  vary 
with  changes  in  chamber  pressure,  thrust  and  nozzle  geometry.  The  propellant 
and  other  motor  parameters  remain  unchanged. 


RX I RL  POSITION  (h) 


Chamber  Pressure:  590  psi;  Thrust:  1000  lbf;  Nozzle  Throat 
Diameter:  10.78  in;  Nozzle  Exit  Diameter:  2.325  in. 

Figure  37.  In-Band  (2. 4-2. 7  urn)  Radiation  Axial  Profiles 
for  15-lb  BATES  Motor  Plumes  with  Identical  Propellants 

In  light  of  this  only  a  qualitative  assessment  of  the  30  firing  cases 
with  IR  emission  data  can  be  done.  Figure  38  depicts  the  typical  plume 
emission  data  at  the  water  band  taken  as  a  function  of  time.  The  chamber 
pressure  curve  is  also  shown  for  reference.  Note  that  at  ignition  the 
radiance  level  pulses  sharply,  as  it  does  near  the  firing  termination.  This 
can  be  attributed  to  the  BKNO3  ignition  pellets  and  initial  flash  at  ignition 
and  the  characteristic  high  pressure  bump  near  burnout.  All  the  data  show 
these  attributes. 

Because  of  its  uncalibrated  form,  not  much  significance  can  be 
attached  to  this  data  beyond  the  presence  of  the  aforementioned  radiance 
spikes.  Note  the  lengthy  H20-band  steady  state  radiance.  This  indicates, 
albeit  after  the  fact,  that  a  stringent  time  resolution  requirement  is  not 
required  in  a  spectral  scan  if  one  assumes  the  same  steady  state  conditions 
exist  in  other  bands  of  interest.  This  fact  is  used  in  planning  the  scanning 
rate  of  future  AFRPL  spectroradiometric  efforts. 
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Figure  38.  Typical  13-lb  BATES  Motor  Water  Band  IR  Emission  Trace 


3.4  Fluae  UV  Signatures 

3.4.1  Data  Reduction  Methods  -  The  fundamentals  of  plume  UV  emission  are  not 

well  known.  With  many  composite  solid  rocket  motor  propellants,  the  gas  phase 

chemistry  is  not  wholly  characterized  or  understood.  There  is  not  much  data 

in  support  of  chemical  modeling  efforts  although  some  spectroscopic  laboratory 

measurements  are  underway  (Ref.  17).  Because  of  these  major  uncertainties  in 

identifying  UV  radiation  mechanisms  for  solid  rocket  propellant  plumes  and 

without  the  imediate  prospect  of  a  supporting  model  to  base  measurements  on, 
o 

2750-4250  A  emission  spectra  at  the  plume  centerline  and  near  the  motor  exit 
plane  were  recorded  for  possible  use  in  a  data  base.  Subsequent  scans  were 
displayed  simultaneously  to  show  temporal  spectral  behavior. 

Reactive  radical  plume  species  expected  to  be  present  in  this  spectral 
region  included  OH,  CN,  NH,  CH,  NO  and  CO.  Some  stability  additives  were  also 
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expected  to  possibly  contribute.  It  was  also  hoped  a  time  history  of  plume 
species  would  provide  useful  information.  Accurate  detection  of  the  species 
present  within  the  plumes  is  highly  dependent  on  the  spectral  resolution  and 
sensitivity  of  the  measuring  instrument,  and  as  will  be  illustrated,  these 
were  constraining  aspects  of  the  AFRPL  UV  measurements. 

Section  2.4.1  describes  how  potentiometer  output  signals  proportional 

to  the  wavelength  position  of  the  UV  monochromator  grating  drive  were  curve- 

fitted  via  the  linear  regression-derived  relation  of  Equation  15.  The 

photomultiplier  voltage  response  to  the  throughputted  UV  signal  was  then 

plotted  as  a  function  of  the  wavelength  for  each  scan  direction.  Computer 

software  was  created  to  plot  successive  scan  directions  as  a  function  of 

wavelength,  yielding  time  contours  of  the  UV  signature.  The  uncertainty  of 

o 

the  wavelength  calibration,  although  different  for  each  firing,  is  about  15  A. 
o 

The  50  A  spectral  resolution  of  the  monochromator  was  the  limiting  precision 
for  wavelength. 

As  with  the  IR-E/A  measurements  reported  in  Section  3.3,  the  lack  of 
an  accurate  absolute  intensity  calibration  (due  to  poor  documentation  and 
uncharacterized  window  and  lens  spectral  transmittances)  prevents  disclosure 
of  radiation  units;  rather,  intensity  heights  are  displayed  in  their 
uncalibrated  voltage  units.  This  leads  to  a  qualitative  interpretation  only 
of  the  magnitude  of  the  UV  radiation  detected. 

* 

3.4.2  UV  Results  and  Interpretation  -  There  were  twenty-six  15- lb  BATES 

...  ° 
firings  with  successfully  recorded  UV  scans  between  2750  and  4100  A.  Three 

reduced  smoke  composite  propellant  series  were  evaluated:  KKN  (AFRPL  240B 

series),  CWB  (AFRPL  241B  series)  and  DWP  (AFRPL  242B  series).  The  propellent 

compositions  are  shown  in  Table  2  (Section  2.2) 

Representative  UV  spectra  from  each  of  these  series  are  shown  in 
Figure  39.  There  is  similarity  between  individual  UV  signatures  within  each 
series,  but  there  are  significant  differences  between  the  propellant  series  UV 
spectra,  as  illustrated  in  Figure  40. 

Because  of  the  variation  of  several  ingredient  loadings  between  the 
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Figure  39.  Voltage  Re8ponae  va  Wavelength  for  16  Mid-Firing  Half-Scans 
for  KKN  Propellant  (Firing  240B-004) 

three  propellants,  it  is  difficult  to  assert  a  quantitative  cause  for  the 
signature  differences,  but  some  general  observations  can  be  made.  First  and 
most  obviously,  the  time  variation  of  relative  UV  signatures  can  be  gauged 
from  Figures  39-41.  The  integrated  scan  signatures  can  be  seen  to  change 
slightly  through  the  firings.  Some  changes  during  the  steady  state  firing  are 
to  be  expected  because  of  spurious  noise  signals  mainly  due  to  small  plume 
radiance  fluctuations.  As  shown  in  the  UV  spectra,  most  of  the  scan-to-scan 
variations  are  minor  and  can  be  ascribed  to  these  fluctuations.  The  only 
significant  change  between  half-scan  signals  is  shown  in  Figure  39.  This 
sudden  post-steady  state  intensity  growth  was  seen  to  occur  periodically  in 
all  the  three  propellant  series.  A  UV  signal  "burst"  occasionally  appears 
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Figure  40.  Voltage  Response  vs  Wavelength  for  16  Mid-Firing 
Half-Scans  for  CWB  Propellant  (Firing  241B-004) 


near  the  end  of  the  burn  time.  It  is  difficult  to  explain  this  transient 
behavior;  it  may  arise  because  of  an  uneven  mass  flux  of  UV  radiating  species. 
Notice  this  is  a  continuum  effect.  Similar  intensity  versus  time  fluctuations 
have  been  recorded  before  (Ref.  20). 


Figure  42  makes  especially  evident  the  presence  of  some  type  of 
continuum  radiation.  In  this  figure  there  is  a  general  gradual  increase  in 
the  response  of  the  detector  with  increasing  wavelength.  Atop  this  pedestal 
are  signal  fluctuations  characterizing  the  molecular  band  structure  of  plume 
specie.  This  continuum  base  has  been  subject  to  previous  analysis  (Ref.  21). 
Two  mechanisms  for  its  appearance  have  been  postulated:  thermal  greybody 

emission  from  hot  alumina  particulates  and  chemiluminescence,  primarily  from 
the  CO  +  0  reaction,  from  plume  afterburning.  The  relative  contributions 
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Figure  41.  Voltage  Response  vs  Wavelength  for  16  Mid-Firing 
Half-Scans  for  DWP  Propellant  (Firing  242B-004) 

from  these  phenomena  is  uncertain,  but  it  appears  that  both  are  of  some 
importance  in  explaining  the  observed  rocket  plume  UV  emission. 

The  interesting  difference  the  data  presented  in  this  report  poses  is 
that  none  of  the  three  propellant  formulations  shown  in  Figures  39-42,  KKN, 
CWB  and  DWP  had  particulate  additives;  this  is  shown  in  Table  2  (Section  2.2). 
It  is  suspected  then  that  the  observed  continuum  behavior  arose  mainly  from  CO 
+  0  chemiluminescence  rather  than  thermal  emission.  It  is  unfortunate  that 
the  UV  monochromator  was  not  completely  and  reliably  calibrated,  as  the  extent 
of  the  apparent  chemiluminescence  would  be  of  interest.  There  is  a  caution  to 
the  chemiluminescence  assumption,  however:  although  the  UV  monochromator's 

FOV  was  carefully  aligned  to  exclude  the  motor  nozzle,  thermal  emission 
emanating  from  this  source  could  be  a  contributory,  albeit  doubtful,  continuum 
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Figure  42.  Single  Half-Scan  Comparisons  of  Representative  Voltage 
Responses  vs  Wavelengths  for  Three  Reduced  Smoke  Propellants 


radiation  factor.  Other  radiation  sources  may  be  unaccounted  particles  within 
the  propellant  formulation,  out-of-field  scattering  and  carbon  ablation  from 
the  nozzle. 

Neither  the  expected  (1,0),  t2,l)  and  (0,0)  OH  vibrational  bands  at 
o  o 

about  2830,  2890  and  3060  A  nor  the  several  less  than  3000  A  NO  y  bands  were 
found  in  the  data  taken  with  these  BATES  firings.  A  similarly  poor  matchup 
was  found  in  the  UV  data  analysis  contained  in  Reference  32,  although  OH  and 
NO  band  structures  were  reported  with  high  resolution  on  composite  propellant 
strands  in  a  laboratory  environment  (Ref.  21).  One  reason  for  the  poor  AFRPL 
field  measurement  results  was  that  the  sensitivity  of  the  UV  instrument  is 
relatively  low  at  the  correspondingly  low  UV  signatures  of  these  wavelengths. 
The  above  mentioned  laboratory  results  show  that  these  species  radiate  in 
these  regions. 

There  appears  to  be  a  3360  a  NH  emission  peak,  especially  in  the 


latter  scans  of  a  KKN  propellant  firing  (Fig.  39).  The  other  apparent  species 

o 

present  are  violet  CN  bands  which  extend  from  approximately  3500  to  4600  A. 

The  scan  resolution  is  unfortunately  not  fine  enough  to  discern  individual  CN 

.  .  o 

line  features  as  the  3850,  3871  and  3883  A  triplet,  nor  is  it  fine  enough  to 

o 

distinguish  CN  emission  from  an  Fe  molecular  band  at  3860  A.  Pb  and  PbO 

molecular  line  features  characteristically  dominate  the  spectra  from  about 
o 

4060  and  4270  A,  and  this  is  apparently  true  with  the  cases  presented  by  these 

three  propellants.  As  there  is  no  other  explanation,  we  must  assume  trace 

elements  of  Pb  or  PbO  from  the  motor  case  are  in  the  reduced  smoke  AP 

propellant  plumes  tested.  One  final  note:  there  is  an  unidentified  spectral 

o 

triplet  feature  centered  at  about  3760  A  which  is  easily  discerned  with  the 
CWB  propellant  of  Figure  40.  In  the  absence  of  any  other  information,  we  must 
assume  this  is  Fe;  that  is,  in  addition  to  Pb  or  PbO,  trace  elements  of  Fe  are 
also  found  within  the  exhaust  plume  measured. 

4.0  CONCLUSION 


4.1  Summary  and  Discussion 

In  summary,  the  experimental  apparatus  described  herein  was  applied  to 
four  specific  plume  measurement  areas:  visibility  figure  of  merit,  particle 
sizes,  temperatures  and  partial  pressures  and  UV  signatures.  In  each  of  these 
areas,  successful  measurement  is  primarily  contingent  on  instrumentation, 
retrieval  algorithms  and  physical  assumptions. 

The  plume  visibility  figure  of  merit  criterion  used  in  this  study, 
laser  transmissivity,  is  but  one  possible  description.  Drawbacks  such  as  the 
presence  of  obscuring  ignition  smoke,  and  the  monochromaticity  of  the 
measurement  were  enumerated.  The  major  drawback  is  the  lack  of  correlating 
data  between  the  attenuation  measurement  and  total  plume  detectability. 
Nonetheless,  applying  the  transmissivity  figure  of  merit  to  a  series  of 
reduced  smoke  propellant  plumes  indicates  a  difference  between  different 
formulation  exhausts.  By  this  criterion,  a  relative  ranking  of  signatures  can 
be  made.  It  must  be  cautioned,  however,  that  further  work  must  be  done  in 
making  different  figure  of  merit  measurements  and  positively  relating  any 
figure  of  merit  measurements  to  total  plume  detectivity.  Larger  sample  sizes 
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should  also  be  acquired,  both  to  assure  statistical  confidence  intervals  and 
derive  further  figures  of  merit. 

AFRPL  particle  size  measurements  are  dependent  on  a  number  of  physical 
assumptions  and  are  subject  to  difficult  numerical  data  retrieval  methods. 
The  technique  finally  decided  on  in  this  work  yields  an  effective  volume-to- 
surface  mean  particle  diameter  from  the  laser  transmissometer  experiment. 
Results  thus  far  generated  come  from  a  relatively  small  sample  of  solid  rocket 
motor  firings,  so  generating  a  conclusive  relationship  between  the  retrieved 
particle  diameter  and  motor  geometry  or  performance  cannot  yet  be  done. 
Comparing  the  results  with  a  semi-empirical  curve  fit  (based  on  earlier  data 
on  captured  particles)  of  expected  mean  diameters  shows  the  retrieved  particle 
diameters  are  significantly  larger,  but  again,  more  data  must  be  taken  before 
a  modification  to  the  curve  fit  can  be  considered.  Further  work  must  also  be 
done  to  verify  and  calibrate  the  retrieved  particle  diameters  and  to  generate 
an  experiment  and  analytical  technique  to  confidently  deconvolve  particle  size 
distribution  functions. 

Recognized  biases  introduced  by  particle  capturing  techniques  led  to 
the  design  of  the  supersonic  shock-swallowing  particle  collection  probe 
described  within  the  report.  Of  the  collected  particles,  an  interesting 
phenomena  arises:  apparently  hollow  particles  embellished  with  a  hole. 
Higher  effective  emissivities,  particle  nonsphericity  and  a  larger  number  of 
particles  for  a  constant  mass  loading  are  possible  results  of  this. 

Plume  temperatures  and  partial  pressure  proved  analytically  impossible 
to  derive  from  multispectral  scan  results  due  mainly  to  an  insufficient 
emission/absorption  path  variation  with  wavelength  in  particle-loaded  plumes. 
The  scanning  monochromator  was  then  locked  at  the  2.7  pm  H2O  p  branch  band  to 
gauge  the  time  trace  of  H2O  emission  and  absorption.  It  is  clearly  shown  to 
follow  the  pressure  trace,  remaining  practically  invariant  during  the  firing's 
steady  state  portion.  This  information  was  used  in  the  time  resolution 
specification  of  an  AFRPL  procurement  of  a  later  infrared  signature 
measurement  device,  a  Fourier  Transform  Spectrometer.  Unfortunately,  because 
of  the  uncertainty  in  calibration  procedures,  single-wavelength  and  line-of- 
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sight  temperature  and  partial  pressure  retrieval  could  not  be  done  from  the 
emission-absorption  measurement. 

o 

Plume  UV  emission  spectroscopy  was  accomplished  from  2750  to  4100  A 
for  three  reduced  smoke  composite  propellant  series.  The  UV  spectra  time 
evolution  was  shown  with  single  scan  contours  plotted  along  a  time  axis.  An 
observed  continuum  effect  was  due  more  to  chemiluminescence  than  thermal 
emission  because  exhaust  particulates  were  not  present,  with  the  exception  of 
carbon  ablation  from  nozzle  erosion.  Several  expected  OH  and  NO  bands  are  not 
present  in  the  data;  although  NH  and  violet  CN  emission  bands  features  were 
identified.  Trace  Pb  or  PbO  and  Fe  elements  were  perceived  to  be  present 
within  the  UV  signature.  Because  UV  modeling  is  still  in  its  early 
developmental  stages  future  and  better-resolved  UV  spectra  in  similar  field 
applications  would  yield  data  to  supplement  these  results. 

4.2  Future  Work 

The  follow-on  in-house  effort  is  tailored  more  specifically  after 
better  characterizing  the  sensitive  input  files  to  various  JANNAF  signature 
prediction  codes.  The  particle  sizing  effort  will  still  be  pursued  by 
supplementing  the  laser  transmissometer  experiment  with  the  supersonic  shock 
swallowing  quench  probe  described  in  Section  2.3.  Analytical  and  experimental 
particle  sizing  methods  will  continue  to  be  investigated  to  eventually 
retrieve  ensemble  particle  size  distributions  within  the  plume  flowfield.  Of 
more  driving  concern,  however,  are  other  particle  properties,  specifically, 
the  infrared  optical  properties  of  particle  indices  of  refraction  and 
emissivities.  This  dictated  the  supersonic  inlet  probe  design. 

AFAPL  is  in  the  process  of  acquiring  more  field  instrumentation  for 
code  verification.  Instruments  include  a  high  resolution  infrared  Fourier 
Transform  Spectrometer,  an  ultraviolet-visible  flat  field  spectrograph,  plus  a 
host  of  other  infrared  and  ultraviolet  radiometric  and  support  devices.  These 
instruments  will  be  transported  to  different  test  stands  aboard  the  just- 
renovated  Mobile  Radiometry  Trailer  (MRT).  Appendix  D  describes  plume 
measurements  to  be  done  on  15-lb  BATES  motors  before  these  plans  and  MRT  are 
implemented. 
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APPENDIX  A 


A  BRIEF  REVIEW  OF  THE  THEORETICAL  BASIS  OF  THE  AFRPL  ROCKET  EXHAUST  PARTICLE 
SIZING  EFFORT:  ASSUMPTIONS  AND  MODELING 


This  appendix  gives  the  basis  and  assumptions  germaine  to  the  AFRPL 
particle  size  data  reduction  efforts.  Only  the  Mie  scattering  theory  tools 
needed  to  understand  why  the  retrieval  methods  given  in  the  text  of  this 
report  (Section  3.1)  were  used  are  covered  here.  Whenever  concepts  or 
equations  are  given  in  this  appendix,  appropriate  references  with  the 
necessary  detailed  explanation  or  derivations  are  cited.  It  must  be 
remembered  that  the  purpose  of  this  appendix  is  to  present  a  primer  on  the 
foundation  of  light  scattering  theory  only  as  it  applies  to  the  AFRPL  plume 
particle  sizing  laser  transmissometer  experiments  and  not  to  present  a  rigid 
dissertation  of  Mie  theory. 

1.0  ROCKET  PLUME  APPLICATION  ASSUMPTIONS 

A  laser  beam  transmitting  through  a  particle-laden  rocket  exhaust  flow 
is  but  a  specialized  application  of  the  more  general  case  of  electromagnetic 
radiation  traversing  any  volume  of  particulates.  Two  extinction  mechanisms 
are  responsible  for  the  portion  of  energy  removed  with  the  passage  of  the 
laser  electromagnetic  radiation  through  a  particulate-loaded  medium:  particle 
absorption  and  particle  scattering.  These  two  mechanisms  are  highly  dependent 
on  properties  of  the  particle  cloud.  The  properties  include  particle  size, 
shape,  mass  loading,  density  and  refractive  index  (m  *  n  +  ik(X)).  An 

additional  influencing  parameter  is  the  wavelength  X  of  incident  radiation. 
To  model  these  and  other  affecting  parameters  is  a  difficult  undertaking,  one 
which  is  severely  constrained  by  assumptions  in  geometry,  size  distribution, 
plume  density  fluctuations  and  degree  of  multiple  scattering — which  are  in 
turn  affected  by  particle  size,  shape,  mass  loading,  etc.  (Refs.  A-l  and  A- 
2). 
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However,  some  simplifying  assumptions  (toward  an  absorption/scattering 
model  and  inversion  scheme)  can  be  applied  to  the  special  case  the  AFRPL  laser 
transmissometer  experiment  presents.  The  polarized,  coherent  and 

monochromatic  radiation  provided  by  laser  radiation  provides  known  values  for 
the  incident  polarization,  electromagnetic  radiation  coherence,  and  wavelength 
parameters.  Assuming  optically  thin  plumes  further  simplifies  the  model  by 
eliminating  consideration  of  secondary  (or  multiple)  scattering  and  density 
fluctuations  found  with  optically  thick  media.  This  seems  a  reasonable  caveat 
given  the  small  particulate  mass  loading  of  low  visibility  plumes  (Ref.  A-3), 
especially  when  coupled  with  the  small  scattering  focal  volume  defined  by  the 
AFRPL  experimental  design.  To  assume  that  the  particle  size  distribution 
function  (PSDF)  n(r)  (given  in  number  of  particles  vs  size)  is  some  continuous 
function  also  seems  reasonable  given  the  current  state  of  knowledge  of  solid 
rocket  motor  plume  particle  properties  (Ref.  A-4). 

A  final  and  crucial  assumption,  spherical  particles,  is  subject  to 

some  debate.  Figure  A-l  is  a  photomicrograph  of  some  AI2O3  particles  captured 

at  AEDC  from  a  0.25-lb  solid  rocket  motor  which  shows  that  the  assumption  of 

rocket  exhaust  particle  sphericity  may  have  some  basis  (Ref.  A-5).  However, 

surface  irregularities  arising  from  outer  shell  deposits,  cracks  on  the 

surface,  and  agglomeration  in  addition  to  even  slight  nonuniform  particle 

radii  of  curvature  will  render  a  particle  nonspherical,  and  a  cloud  of 

polydisperse  nonspherical  particles  illuminated  by  laser  light  may  yield 

fluctuations  in  the  ensuing  scattered  light  coherence  and  intensity.  It  has 

% 

also  been  contended  that  the  random  orientation  in  such  a  cloud  renders  an 
effective  monodisperse  cloud  of  spherical  particles  with  some  characteristic 
surface-equivalent  diameter,  but  this  postulate  is  thus  far  unsubstantiated 
(Ref.  A-6).  What  experiments  do  show  is  that  Mie  theory  apparently  holds  for 
approximately  monodisperse  nonspherical  and  randomly  oriented  particles  within 
the  size  range  we  are  interested  in  (Ref.  A-7). 

Scattering  parameters  from  pure  geometries  as  spheres,  cylinders  (Ref. 
A-8),  arbitrary  axisymmetric  bodies  (Ref.  A-9),  and  arbitrary  smooth  convex 
surfaces  of  variable  curvature  (Ref.  A-10)  have  been  calculated.  Single¬ 
scattering  from  cubic  particles  has  also  been  measured  (Ref.  A-ll).  However, 
no  rigorous  analysis  of  the  effects  of  discontinuities  on  the  particulate 
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(A-l ) 


where 


I0  *  intensity  of  the  incident  light; 
a  *  2ftr/X  ■  dimensionless  particle  size  parameter; 

^  =  laser  wavelength; 

m  *  the  complex  index  of  refraction, 

9  =  the  aspect  angle  between  I0  and  Is;  and 

A 

F  =  the  normalized  four-by-four  scattering  transformation  matrix 
(often  called  a  Meuller  matrix  since  it  is  found  via  Mueller  calculus). 

Four  Stokes  parameters  comprise  IQ,  and  fully  describe  the  incoming 
electromagnetic  radiation  and  its  polarization  states.  Equation  A-l  shows 
that  it  is  the  vector  of  these  Stokes  parameters  that  is  transformed  from 
incident  to  scattered  radiation  by  the  Mueller  matrix  operator  F  (Ref.  A-7). 


The  Stokes  parameters,  I(a  ,  m,  6  )  Q(a,  m,  0),  U(a ,  m,  6 )  and  V(a,  m, 
6),  can  be  further  reduced  to  two  intensity  distribution  functions  described 
by  their  dimensionless  complex  amplitudes  S}  and  S2.  These  amplitudes  are 
derived  in  Reference  A-7,  where 


|Sll 


E  2n  +1 
n-1  n(n  +1) 


bnTn) 


(A-2) 


giving 


I  SX  |  2  -iRe^x)  +  Im(Si)|2. 


(A-3) 


Similarly, 


This  can  also  be  reduced  to  a  relation  like  Equation  A-3  where 


|  S2 1  2  =■  |  Re(S2)  +  Im(S2)|2.  (A-5) 

In  Equations  A-2  through  A-5,  n  =  order  of  electric  (an)  or  magnetic  (bn) 
particle  waves,  where  an  and  bn  are  complex  functions  derived  from  Bessel 
functions  and  are  dependent  on  a  and  m  but  not  on  8.  Conversely,  the  phase 
functions  and  depend  only  on  6  and  use  first  and  second  derivatives  of 
Legendre  polynomials  of  order  n  and  argument  cos  9. 

^n(z)  and  ^n(z)  are  the  Ricatti-Bessel  functions  used  in  determining 

a  and  b  where 
n  n 


1>  <z)  -  (Vrz)**  J  .  u(z) 
n  n  + 


(A-6) 


and 


«n(z) 


(*5pz)  H 


(2) 


+  % 


(A-7) 


(2) 

Jn  +  and  H  n  +  are  Be88®l  functions  of  the  first  and  third  types.  Given 
these  definitions, 


a  8  il£g2  *  ixSll  ~  a^n(y)^n'(x) 

n  ♦«(y)5_(x)  -  “♦_(y)€'6o  (a-8) 

n  n  n  n 


and 


b  =  -  4>n(y)^  1  (°0 

n  ra<Cy)5n(a)  -  (y) 5^(a)  (A-9) 


where,  as  defined  earlier, 


a  a  2irr/X  ; 


m  *  complex  index  of  refraction;  and 


y  ■  mot  . 
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Reference  A- 2  provides  a  derivation  for  Equations  A-8  and  A-9. 


Also  appearing  in  Equations  A-2  and  A-4  are  *n  and  Tn«  Expressed  in 
terms  of  first  and  second  derivatives  of  the  Legendre  polynomials  p,  they  are 


/  _  dP„(cos  0) 

nn<cos  0)  “  ~~d  cos"  9 


(A-10) 


t  (cos  9)  *  cos  *n  (cos  6)  -  sin29 
n  11  a  cos  0 


(A-l 1 ) 


3.0  SINGLE  PARTICLE  EFFICIENCY  FACTORS 


We  can  define  efficiency  factors  for  scattering  Q8C^®»  ■)  and 
extinction  Qe(a,  m).  Efficiency  Q  is  the  ratio  of  energy  removed  by  particle 
extinction  or  scattering  to  its  geometric  cross  section.  The  Qbc  relation  is 
written  as 


sc  2 

irr 


(A-l 2) 


where 


a  *2irf  o(6)sin0d8 
P  Jo  P 


(A-13) 


and  a  represents  the  total  scattering  cross  section.  The  differential 
P 

scattering  coefficient  is  defined  as 


op(e)  -  *rrmaxd_2M.  r2n(r)dr 
J  rmin  “ 


(A-14) 


where  0  is  the  unit  solid  angle  of  scattered  light.  Using  the  a  relation 
given  in  Equation  A-13,  Equation  A-12  becomes 


2  (v 

~2  °p<6 

t  Jo 


)  sin 8  d  0 


(A-15) 


and  can  be  expressed  in  discrete  form  aa 
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Similarly, 


Qe(a  ’  n)  " 


(2n  +  1) 


T:Re  (a  +  b  ) 
n  n 


(A-17) 


Since  energy  ia  conserved,  the  absorption  coefficient  Qai,8^a»m^  can  be 
expressed  as 


Qabg(<*>  =  Qe(a»  “  Q8C(a»  m)‘ 


(A-18) 


From  Equations  A-12  and  A-14,  the  differential  scattering  cross 

section  a  (6)  is  seen  to  govern  the  amount  of  scattered  light  as  a  function  of 
P 

scattering  angle.  This  can  be  made  explicit: 


E(e)  =  a  (6)  I  (W  cm  ^  sr-*). 
P  ° 


(A-19) 


E(6)  is  the  differential  scattered  power  per  steradian  for  each  unit 
scattering  cross  section. 


An  application  of  the  Meuller  transformation  matrix  to  Equations  A-l , 
A-14  and  A-18  yields 


X^I  rr 


E(X,  m,  8) 


o  /  max  Cl  .  . .  ,  .  . 

—  I  F  (a,  m, 8)n (r)dr , 


(A-20) 


where 


E( x,  m,  9)  ■  the  angular  single-scattered  power  from  a  cloud  of 

particles, 


F( x»  m,  9)  ■  a  normalized  scattering  matrix,  and 


n(r)  *  some  mean  value  for  the  PSDF  (see  Ref.  A-7). 


Equation  A-20  is  recognized  as  a  Fredholm  integral  of  the  first  kind. 
Applying  a  bistatic  numerical  inversion  solution  to  solve  for  n(r)  and  to 
account  for  both  X  and  8  is  nontrivial. 
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4.0 


SUMMARY 


This  appendix  reviews  the  tenants  of  Mie  scattering  theory  applicable 
to  the  AFRPL  particle  sizing  experiment  and  numerical  deconvolution 
approaches.  The  theory  behind  three  particle  sizing  techniques  tried  at 
AFRPL — laser  polarization,  angular  intensity  and  transmission/extinction — was 
presented  in  this  appendix,  albeit  without  much  mention  of  corresponding 
numerical  techniques.  This  is  the  intention;  the  numerical  methods  are  given 
in  the  course  of  the  report  (Section  3.1).  Underlying  assumptions  pertaining 
specifically  to  the  AFRPL  Mie  scattering  experiment  were  also  given  to  further 
constrain  the  extent  of  this  review.  Full  descriptions  and  derivations  of 
light  scattering  theory,  while  interesting,  are  beyond  the  intention  of  the 
appendix;  they  are,  however,  included  in  the  references. 
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APPENDIX  B 


A  BRIEF  REVIEW  OF  THE  THEORETICAL  BASIS  OF  THE  AFRPL  INFRARED 
EMISSION/ABSORPTION  EXPERIMENT  TO  RETRIEVE  PLUME 
SPECIES  TEMPERATURES  AND  PARTIAL  PRESSURES: 
ASSUMPTIONS  AND  MODELING 


This  appendix  addresses  Che  underlying  theory  behind  Che  AFRPL  BATES 
plume  infrared  emission/absorption  (IR-E/A)  experiment.  Neither  a  rigid 
derivation  of  applicable  radiation  models,  a  primer  on  the  numerical  inversion 
techniques  nor  a  description  of  the  apparatus  employed  with  the  AFRPL 
experiment  is  partaken  here;  radiation  models  are  derived  and  discussed  in 
detail  elsewhere  (Refs.  B-l  and  B-2),  the  numerical  deconvolution  procedure  is 
given  in  the  text  of  this  report  (Section  3.3)  and  the  instrumentation 
configuration  is  also  described  within  the  report  (Section  2.4).  The  intent 
with  this  appendix  is  instead  to  familiarize  the  reader  on  the  radiative 
transfer  relationships  used  to  determine  values  for  plume  species  temperature 
and  partial  pressure  profiles  through  application  of  emission  and  absorption 
spectroscopy.  Three  progressively  realistic  radiative  transfer  models  will  be 
discussed:  the  particle-free  homogeneous  gas  case,  the  particle-free 
inhomogeneous  axisymmetrc  gas  case  and  the  particle-laden  inhomogeneous 
axisyimnetric  gas  case. 

1.0  TEMPERATURE  AND  PARTIAL  PRESSURES  FROM  A  PARTICLE-FREE  HOMOGENEOUS  GAS 

This  is  the  simplest  case  for  plume  IR-E/A  analysis,  although  the 
restrictions  imposed  by  assuming  the  absence  of  radiating  and  scattering 
particles  and  also  assuming  a  plume  gas  column  with  no  density  fluctuations 
make  the  model  unrealistic  for  particle-laden  flows.  Nonetheless,  the  first- 
order  species-independent  analysis  will  familiarize  the  reader  with  the 
necessary  basis  for  the  more  stressing  band  model  inversion  schemes. 

If  the  plume  is  at  an  equilibrium  temperature  Tp,  its  monochromatic 
radiance  LP(X,  Tp)  is  given  by  the  relation, 


LP(X,TP)  -ePL®B(X,Tp) 


(  B— 1 ) 


where  the  wavelength  X  and  plum  temperature  Tp  are  constant  and  where: 
3  monochromatic  emissivity  of  the  plume  and 


BB  d 

L“  (A  ,  Tp)  3  monochromatic  blackbody  radiance. 


This  blackbody  radiance  is  governed  by  the  Planck  function, 


L“(X  ,  Tp)  -  - 

X  x5(ec2AT'>  -  1) 


(B-2) 


where 


C1 


2irc2h  -  3.7419  X  10-12  Wcm-2  and 


C2  -  hc/Kb  -  1.439  cmK  . 


Furthermore,  Kirchoff'a  law  defines  the  equivalence, 


(B-3) 


where 


aP  3  monochromatic  absorptivity  of  the  plume  and 

A 

xP  3  monochromatic  transmissivity  of  the  plume. 

A 

Combining  Equations  B-l,  B-2  and  B-3  yields  an  expression  for  the 
plume  temperature 


In  °1  (1  -  xP) 
LPX5  -  1 


(B-4) 


as  a  function  of  monochromatic  wavelength  X,  plume  transmittance  tp  and  plus 
radiance  LP.  Note  that  TP  is  independent  of  X. 


According  Co  Che  Beer-Lambert  principle,  a  specific  relaCion  exisCs 
beCween  transmittance  and  partial  pressure  for  purely  absorbent  low 
concentration  media,  i.  e., 

In  TP  “  -cpPpKp(c,  p,  T)L  (B-5) 

A  A 


where 

cp  *  concentration; 

P  *  pressure  of  the  plume; 

Kp(c,  p,  T)  “  monochromatic  absorption  coefficient  of  the  plume  (cm-* 

A 

atm-* ) ;  and 

L  3  length  of  transmitting  path  s  through  any  arbitrary  plume  line  of 
sight. 

The  plume  monochromatic  concentration  cp  is  easily  derivable  from 
Equation  B-5,  given  some  input  value  for  static  pressure.  In  the  least  case, 
the  partial  pressure  can  be  extracted  from  the  above  equation. 

2.0  TEMPERATURE  AND  PARTIAL  PRESSURES  PROFILES  FROM  A  PARTICLE-FREE 
INHOMOGENEOUS  AXISYMMETRIC  GAS 


For  a  line  of  sight  s  through  a  nonuniform  and  gaseous  plume  medium 
containing  temperature  and  pressure  gradients  as  illustrated  in  Figure  B-l, 
the  resulting  monochromatic  plume  radiance  Lp  can  be  expressed  in  simplified 
form  as 


Lp( X,  Tp) 


-  f  (LBB(s)  +  LE(s)) 
-0 


drP (s) 
d  s 


ds 


(B-6) 


where 


L  3  distance  along  line  of  sight  s, 

BB 

L  (s)  *  monochromatic  blackbody  radiance  at  s, 

A 
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HOT  GAS  REGION 


t 


Figure  B-l.  Gas-Only  Flow  Optical  Geometry 


[f 

L  (s)  •  monochromatic  external  radiance  incident  on  plume  boundary  and 


t*(s)  *  average  monochromatic  transmittance  of  the  plume  at  s. 


Substituting  Equation  B-2  into  Equation  B-6  and  solving  for  the  plume 
temperature  profile  Tp(s)  requires  some  numerical  inversion  to  solve  for 


(s)  (Ref.  B-2).  Tp(s)  is  again  independent  of  *. 


Using  the  Beer-Lambert  absorption  law,  transmittance  can  now  be 
expressed  as 


'  ■  ; 

InxP(s)  ■/  cp(s)pP(s)Kp(s)y^(s)ds. 


(B-7) 


Kp( s)  *  monochromatic  absorption  coefficient  of  the  plume  and 


yj(s)  *  band  model  function. 


Equation  B-7  differs  from  Equation  B-5,  as  terms  are  written  as 
functions  of  the  line  of  sight  s;  this  is  a  concession  to  the  plume 
inhomogeneity. 


The  band  model  function  yp(s)  is  necessarily  inserted  in  Equation  B-7 


to  describe  plume  species  band  system  structure  since  they  are  highly 
sensitive  to  fluctuations  in  temperature,  partial  pressure  and  gas  density. 


 V.  k  *.  V,  Vi  v*  Vi  *  .  k  Vs  V.s  4  •  v 1  L  «ikA  • 
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The  band  model  must  have  line-by-line  resolution  to  simulate  the  plume 
species'  radiative  transfer  (most  notably  with  the  primary  combustion  gases, 
H2O  and  CO2)  because  of  these  sensitivities.  Available  band  models  fit 
empirical  data  very  well  (Refs.  B-3  and  B-4).  The  particular  function  y£(s) 
used  here  arises  from  the  Malkmus  statistical  band  model  and  the  calculations 
are  carried  out  with  a  Curtis-Goodson  nonuniformity  approach  fully  described 
in  Reference  B-4.  This  function  accounts  for  the  difference  between  the  Lund 
and  Beer-Lambert  absorption  models  and  is  highly  dependent  on  postulated 
spectral  line  shapes. 

The  numerical  method  to  retrieve  concentration  and  pressure  profiles 
from  a  model  reflected  by  Equation  B-7  is  discussed  in  the  text  of  this  report 
(Section  2.4)  and  detailed  in  Reference  B-2. 

3.0  TEMPERATURE  AND  PARTIAL  PRESSURE  PROFILES  FROM  A  PARTICLE-LADEN 
INHOMOGENEOUS  AXISYMMETRIC  GAS 


Coupling  the  gas-only  observed  radiance  relation  (Equation  B-6)  with 
single  scattering  effects  of  particles  within  the  plume  flowfield  yields  a 
revised  monochromatic  plume  and  radiance  LP  expressed  as 


.L  _  r  4tt 

LP(X,TP)  =  j  (LBB(s)  +  LE(s))  -d'T^5--  dsj  a 
_  n  J  n 


(s)  S(s)d0 


(B-8) 


where 


a(s)  ■  particle  absorption  cross  section  and 


g(s)  *  particle  scattering  cross  section. 


These  cross  section  terms  are  integrated  over  the  solid  angle  defined  by  the 
scattering  angles  6  and  $  shown  in  Figure  B-2. 


Equation  B-7  can  be  likewise  modified  to  include  particle  effects. 
The  relation  governing  plume  transmittance  for  a  two-phase  system  is 


lnxp(s)  •‘•J  cp(s)pp(s)kp(s)y£(s)ds  J  a(9)&(s)ds. 


(B-9) 


The  numerical  algorithm  to  be  used  with  Equation  B-9  ia  very  similar 
to  that  for  Equation  B-7.  Theory  and  premises  behind  the  single  particle 
light  scattering  theory  used  in  Equations  B-8  and  B-9  are  discussed  in 
Appendix  A  of  this  report. 


^  /A 

//'&' 
_ y  v 


SENSOR 


Figure  B-2.  Gas-Particle  Flow  Optical  Geometry 
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APPENDIX  C 


A  BRIEF  REVIEW  OF  THE  ADAPTIVE  NOISE  CANCELLING  SIGNAL 
PROCESSING  TECHNIQUE  AS  APPLIED  AT  AFRPL 

Difficulty  in  reducing  the  noisy  scattering  detector  data  was  cited  in 
the  text  of  the  report  (Section  3.1)  as  a  practical  experimental  problem 
encountered  in  addition  to  the  analytical  and  computational  difficulties 
encountered  with  the  particle  size  distribution  retrieval  algorithm.  Despite 
the  deleterious  effects  of  not  having  an  accurate  size  distribution  retrieval 
code,  a  signal  processing  scheme  to  overcome  low  scattering  detector  signal- 
to-noise  ratios  was  still  pursued.  A  major  reason  for  this  was  to  provide  a 
potentially  useful  processing  technique  for  future  AFRPL  data  reduction 
efforts  (e.  g.,  as  described  in  the  report  under  Section  3.2). 

This  appendix  describes  the  concepts  and  application  of  a  signal 
processing  technique  called  adaptive  noise  canceling  and  shows  result  of  its 
application  to  data  from  a  scattering  detector  the  signal  test  run. 

1.0  THE  CONCEPT  OF  ADAPTIVE  NOISE  CANCELLING  (REF.  C-l) 

Figure  C-l  describes  the  structure  of  the  adaptive  noise  canceling 
system  used  in  the  AFRPL  application  to  reduce  the  additive  noise  that 
corrupts  desired  scattering  detector  signals.  The  noise  cancellation  can  be 
carried  out  in  real  time  or  during  the  post  processing  of  the  recorded  signal. 
For  the  system  to  work  effectively  the  following  conditions  are  necessary: 


Figure  C-l.  An  Adaptive  Noise  Cancelling  System 


(1)  the  desired  signal  s  is  corrupted  by  an  additive  noise  n^ 
uncorrelated  with  the  signal  s;  and 

(2)  availability  of  a  noise  n^  uncorrelated  with  s  but  correlated  in 
some  unknown  way  with  the  noise  Oq. 

Note  that  for  the  signals  encountered  with  the  scattering  detector  case,  (1) 
is  a  reasonable  assumption,  and  the  noise  signal  n^  of  condition  (2)  can  be 
recorded  during  the  firing. 

The  idea  of  the  adaptive  noise  canceler  is  to  filter  the  noise  signal 
n^  to  produce  an  output  y  that  is  as  close  a  replica  as  possible  of  n^  by 
dynamically  and  automatically  adjusting  the  filter's  own  impulse  response. 
The  adjustments  are  carried  out  through  an  algorithm  that  responds  to  an  error 
signal  dependent,  among  other  things,  on  the  filter's  output.  Thus,  with  a 
proper  algorithm  the  filter  can  operate  under  changing  conditions  and  can 
readjust  itself  continuously  to  minimize  the  error  signal  in  the  Least  Mean 
Square  (LMS)  sense. 

To  justify  how  the  system  works  assume  that  s,  nQ,  n^  are 
statistically  stationary  and  each  has  a  mean  of  zero.  Consider  the  system 
output  z: 


z  -  s  +  nQ  -  y.  (C-l) 

Squaring  and  taking  expectations  gives: 

E(z2)  -  E(s2)  +  E(nQ  -  y)2  ♦  2E(s(n0  -  y)).  (C-2a) 

But  E(S(ng  -  y))  =*  0  because  s  is  uncorrelated  with  n^  and  with  y.  Hence, 

E(z2)  »  E(s2)  +  E(n0  -  y)2.  (C-2b) 

2 

From  Equation  (C-2),  note  that  the  independent  signal  power  E(s  )  will  be 

2 

unaffected  as  the  filter  is  adjusted  to  minimise  the  output  power  E(z  ). 
Accordingly,  the  minimum  output  power  can  be  expressed  as 
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(C-3) 


min  (E/z^))“  E(s^)  +  min  (ECn^  -  y)^). 

2  2 

Observe  that  when  the  filter  is  adjusted  so  that  E(z  )  is  minimized,  E(nQ  -y) 
is  also  minimized.  The  filter  output  y  can  be  then  considered  a  best  least 
square  estimate  of  the  primary  noise  nQ. 

The  output  z  will  contain  the  signal  s  plus  noise.  From  Equation  C-l, 

2 

the  output  noise  is  given  by  nQ  -  y.  Since  minimizing  the  output  power  E(z  ) 
also  minimizes  the  output  noise  power  E(oq  -  y)  ,  minimizing  the  total  output 
power  minimizes  the  output  noise  power.  Since  the  signal  in  the  output 
remains  constant,  minimizing  the  total  ouput  power  therefore  maximizes  the 
output  signal-to-noise  ratio. 

2.0  APPLICATION 

To  determine  the  usefulness  of  the  adaptive  noise  canceling  concept 
for  AFRPL  scattering  detector  data,  we  processed  data  from  a  scattering 
detector  channel  (SD1)  taken  in  BATES  242B-015  firing. 

Figure  C-2  shows  the  signal  trace  from  SD1  collected  with  AFRPL  data 
Acquisition  and  Control  System  (DACS).  The  trace  clearly  shows  a  dominant 
effect  of  the  noise  in  the  post-ignition  data. 

To  implement  the  adaptive  scheme  shown  in  Figure  C-l,  a  reference 
signal  correlated  to  the  noise  in  the  SD1  signal  should  be  recorded  during  the 
entire  test.  But  in  the  AFRPL  laser  scattering  experiment  such  signals  were 
not  recorded  because  applying  the  scheme  to  the  data  was  not  considered  at  the 
time  of  the  experiment.  Under  the  situation,  we  generated  an  a-posteriori 
reference  noise  signal  based  on  the  following  assumptions: 

(1)  the  pre-ignition  signal  shown  in  Figure  C-2  is  a  noise  that 
corrupts  additively  the  desired  signal;  and 

(2)  the  pre-ignition  noise  characteristic  remains  the  same  during  the 
entire  firing. 


Using  the  above  assumptions*  we  generated  the  reference  noise  signal 
shown  in  Figure  C-3  by  a  periodic  extension  of  the  pre-ignition  signal  which 
was  long  enough  to  match  the  entire  SDl  signal. 

Figure  C-4  shows  the  adaptive  noise  canceler  output.  As  can  be  seen, 
the  adaptive  filter  rapidly  adapts  to  the  noise  to  produce  zero-level  pre¬ 

ignition  signal  and  continues  working  in  the  post-ignition  period  to  produce 
output  signal  that  is  considerably  different  from  the  original  SDl  signal. 

It  is  impossible  to  judge  quantitatively  how  good  the  result  is 
because  we  do  not  know  for  certain  how  realistic  the  a-posteriori  reference 
signal  is.  However,  if  in  the  future  we  plan  to  use  such  a  signal  processing 
technique  from  the  start  of  an  experiment,  we  can  record  a  more  realistic 

noise  signal  correlated  to  the  corrupting  noise  and  obtain  a  final  signal  with 
much  less  noise.  Although  we  have  only  shown  results  from  one  scatter 

detector  channels  there  is  a  dramatic  signal-to-noise  improvement  over  the 
unconditioned  signal. 
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APPENDIX  D 


DESCRIPTION  OF  THE  AFRPL  PLUME  RADIOMETRIC  MEASUREMENTS 

PLANNED  FOR  1984 


This  appendix  addresses  the  few  follow-on  and  ancillary  plume 
radiometric  measurements  to  be  done  at  AFRPL  before  the  major  follow-on 
project  described  in  Section  4.2.  Several  reduced  smoke  GAP  polymer 

propellants  are  scheduled  to  be  fired  at  the  same  Pad  5a  of  Experimental  Area 
1-32  where  the  other  measurements  were  taken.  A  radiometer  has  been  added  to 
the  AFRPL  plume  instrumentation  ensemble  described  in  the  report  to  measure 
plume  afterburning.  Once  these  measurements  are  taken  they  are  to  serve  as 
the  basis  for  an  AFRPL  plume  infrared  signature  data  base  to  be  used 
eventually  to  evaluate  plume  infrared  signature.  Follow-on  efforts  mentioned 
in  Section  2.2  of  the  report  will  use  more  instruments  to  measure  plume 
infrared  signatures,  providing  a  nyumber  of  options  including  a  variety  of 
aspect  angles,  field  of  views  (FOVs),  optical  bandpasses  and  spectral 
resolutions. 

For  this  particular  effort,  however,  a  Barnes  Model  12-880  infrared 
radiometer  filtered  at  a  2.395  -  2.715  ym  passband  but  with  a  full  width-half 
maximum  (FWHM)  bandpass  of  2.490  -  2.665  ym.  This  FWHM  bandpass  falls  within 
the  nonisobetic  H2O  spectrum  between  2.38  -  2.63  ym  identified  by  Young  (Ref. 
D-l)  for  the  elevated  temperatures  and  pressures  encountered  with  BATES-type 
motors.  Hence,  the  radiometer  is  essentially  a  plume  constituent  H2O  sensing 
device. 


The  radiometer  has  a  pyroelectric  detector  with  a  chopper  and  internal 
reference  cavity  design.  The  10-90Z  response  time  of  this  system  is  50  ms. 
The  normalized  spectral  response  with  the  detector  and  2.395  -  2.715  ym  filter 
combination  is  shown  in  Figure  D-l.  These  are  the  manufacturers'  data. 

The  reason  for  pointing  the  radiometer  at  the  plume  afterburning 
region  is  that  most  of  the  plume's  radiation  emanates  from  this  area,  and  the 
radiometer's  2°  FOV  combined  with  the  limited  space  within  the  test  pad 
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necessarily  constrained  the  measurement  to  only  a  certain  portion  of  the 
plume.  The  radiometer  was  positioned  20.2  feet  from  the  BATES  motor  exit 
plane  and  oriented  at  a  42°  aspect  angle  from  the  motor  centerline,  away  from 
the  plume  flow.  The  afterburning  region  was  identified  by  applying  the  JANNAF 
Standardized  Infrared  Radiation  Model  (SIRRM)  plume  predictive  code  to 
individual  BATES  motor  propellant,  nozzle  and  flowfield  inputs,  and  running 
the  filter  integration  bandpass  module  in  the  spatial  radiant  intensity  mode 
(Ref.  D-2) .  Figure  D-2  is  an  example  of  the  SIRRM  code  output  when  run  this 


Pointing  the  radiometer  is  a  relatively  simple  operation,  done  with  a 
high  intensity  source  placed  at  the  plume  centerline.  In  the  case  with  Figure 
D-2,  the  source  would  be  placed  2  m  downstream  from  the  motor  exit  plane. 


Figure  D-2.  In-Band  Radiant  Intensity  vs  Axial  Position  Along 
15-lb  BATES  Motor  Plume  Centerline 


At  this  writing  no  reportable  data  had  been  taken  via  the  radiometer 
channel  for  the  five  reduced  smoke  static  firings  due  to  problems  in  the 
radiometer's  output  electronics  and  data  acquisition  system.  However,  at 
least  15  more  test  firings  plus  other  AFRPL  targets  of  opportunity  are 
expected  to  occur  soon.  The  following  calibration  section  addresses  how  pre¬ 
firing  preparations  are  made. 


CALIBRATION 

The  calibration  techniques  applied  to  these  measurements  are  detailed 
in  Reference  D-3.  Essentially,  an  extended  blackbody  source  calibration 
technique  was  used  to  calculate  the  calibration  transfer  function  in  the 
absolute  responsivity  calibration.  Field  of  view  mapping  and  background  and 
noise  level  determination  steps  were  undertaken  in  addition  to  consideration 
of  out-of-bound  rejection.  These  error  sources  plus  those  enumerated  in 
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Section  2.4.3  of  this  report  are  embedded  in  a  calibration  uncertainty  factor. 
This  uncertainty  factor  gives  the  total  bias  error  bar  associated  with  the 
radiometer  data.  Because  each  of  the  firings  was  unique,  no  random  error 
uncertainties  could  be  attached  to  the  total  error  margin. 
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